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Observations  of  positively  charged  particles  that  are 


generated  on  or  near  satellite  surfaces  have  been  made  on 


several  spacecraft.  This  thesis  postulates  sputtering  of 


the  satellite  surface  due  to  ambient  ion  impact  as  the 


generating  mechanism.  Calculations  are  made  using  the 


S igmund-Thompson  sputtering  theory  to  determine  the  response 


at  the  particle  detectors.  These  calculations  indicate  that 


surface  sputtering  creates  a  sufficient  flux  to  account  for 


the  observed  phenomena.  The  NASA  Charging  Analyzer  Program 


was  run  to  determine  the  trajectories  for  actually  observed 


particles.  The  calculated  trajectories  were  determined  to 


lead  to  the  spacecraft  surface,  again  indicating  that 


surface  emission  was  the  source.  The  sputtering  flux  as 


calculated  was  insufficient  to  cause  any  significant  short¬ 


term  damage  to  the  spacecraft,  beyond  thin  coating  erosion. 
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A.  OVERVIEW 

Spacecraft  charging  is  the  development  of  a  potential 
difference  between  a  space  vehicle  and  its  plasma 
environment.  It  is  of  significance  for  scientific  pcrpcsss, 
since  the  potential  may  affect  measurement  of  the  ambient 
properties  of  the  environment,  and  for  practical  purposes, 
since  it  may  lead  to  anomalous  command  signals  and  physical 
damage  to  the  affected  spacecraft.  Although  predictions  of 
spacecraft  charging  had  been  made  previously,  the  first 
reports  of  high  level  spacecraft  charging  were  those  for 
Applications  Technology  Satellite  (ATS)  5  {Ref.  1).  Similar 
events  were  subsequently  noted  on  ATS-6  and  P78-2  (CCATHA) 
by  various  observers  (Ref. 2}. 

Generally,  negative  charging  events  are  identified  by  an 
intense  peak  in  the  measured  ion  spectra  at  a  specific 
energy,  with  an  absence  of  ions  of  lesser  energies.  This 


intense  pea) 


most  clearly  seen  on  the  plot  of  the 


distribution  function  versus  energy,  as  will  be  displayed 


later  in  this  work.  This  peak  is  known  as  the 


:harg  inq 


peak.  Guch  observations  result  from  the  acceleration  of  the 
ambient  plasma  ions  into  the  detector  by  the  potential  drop 
between  the  plasma  and  the  spacecraft.  Occasionally,  ion 
fluxes  are  seen  at  energies  below  the  peak  energy,  which 


4 

5? 


would  appear  to  violate  conservation  of  energy  for  a 
CO  1 1  i  s  ’  -  :i  1  a  S3  plasma.  These  ion  d  i  s  tr  i  t  ut  i  ons  have  been 
reported  previously,  and  tentatively  ascribed  to  cputtering 
from  spacecraft  surfaces  {Ref.  3}.  They  have  been  termed 
"spacecraft  generated  ions"  (Ref.  2}. 

Analysis  of  this  particular  phenomenon  is  of  practical 
interest,  as  these  ions  may  contaminate  the  charging  peak, 
leading  to  an  underestimation  of  the  satellite  potential. 
The  contamination  of  the  charging  peak  may  be  particularly 
Important  active  experiments  in  low  earth  orbit, such  as 
SFACELA3.  Additionally,  since  we  believe  the  source  to  t'- 
sputtering,  the  flux  of  these  particles  is  an  indication  .of 
the  damage  to  the  satellite  surface  caused  by  the  vehicle's 
environment.  Also,  the  life  expectancy  of  precision 
satellite  surfaces,  e.g.  optical  surfaces,  and  surface 
coatings  may  be  directly  affected  by  the  sputtering  rate 
appropriate  to  their  environment. 


3.  SCATHA  SPACECRAFT 

Most  of  the  experimental  data  and  calculations  performed 
in  this  thesis  are  based  on  the  P78-2  (SCATHA)  spacecr.aft. 
The  construction,  environment,  and  sensors  of  this  vehic''.  e 
will  be  discussed  in  the  next  sever-al  sections. 

1 .  Spacecraft  Construction 


The  U.o.  Air  Force  P78-2  cp 


W.S  s  1  a  u  n  c  h  e  cl  i 


Tanuary  1379  to  collect  data  for  a  study  of  .hpacecraf 


Ci-iarging  AT  High  Altitudes  (SCATHA), 


in  .3  joint 


he  satellite  body  was 


h'ACA.'Depar tment  of  Defense  Program, 
cylindrical  in  shape  with  a  length  and  diameter  of 

approximately  1 . 7 f  meters.  Seven  booms  were  deployed  in 
orbit  from  the  spacecraft  to  provide  isolation  for 

experiments  from  charging  effects  on  the  satellite  surface. 
Figure  1  depicts  the  SCATHA  spacecraft.  The  top  of  the 
satellite  in  the  figure  is  the  forward  end,  and  the  bottom 
is  aft.  The  experiments  and  commun  i  cat  i  ons  equip.ment  for 
Lhe  satellite  are  primarily  located  on  equipment  decks  in 
the  center  of  the  cylinder.  The  University  of  California, 

San  Diego,  charged  particle  detector  was  on  the  forward  end, 

and  the  ion  gun  was  on  the  aft  end;  both  am  described 
below.  The  outer  cylinder  surface  was  divided  into  three 
general  areas;  two  solar  arrays,  one  forward  and  one  aft, 
and  a  bellyband  between  them  to  provide  access  panels. 
A'^ditionally,  sections  of  the  surface  were  covered  in 
various  materials,  such  as  kevlar  and  mylar,  for  use  in  the 


exper iments 

.  {Ref.  4} 

The 

spacecra  ft 

was 

constructed 

of  spe 

c  ial 

1 ightwe ight 

materials. 

as  i  ts 

orbit  was  at 

the  max 

imum 

capabi 1  ity 

of  the  launch 

vehicle 

The  central 

tube  of 

4-  U  - 

U  i  t 

spacecraft 

was  magnesium  and 

the  equipment 

decks 

were 

s  1  urn i .num  . 

The  solar  a 

rrays 

were  aluminum 

core  wit 

u  - 

fiberglass 

( S i02 )  outer 

face  . 

The  bellyband 

panels 

we  r  e 

covered  with  thermal  paint  and  second  surface  mirrors,  and 


sc  10  2 


acted  as  waste  heat  radiators.  The  forward  cylinder  end  was 
coated  with  gold,  except  for  some  sample  patches.  {Ref.  4} 

2 .  Spacecraft  Orbit  and  Environment 

The  vehicle  was  inserted  into  a  near -geosynchr onous , 
elliptical  orbit  with  a  perigee  of  5.5  R«  and  an  apogee  of 
7.7  R« .  The  orbit  period  was  23.5  hours,  and  the 
inclination  to  the  equatorial  plane  was  7.8  degrees.  The 
satellite  axis  was  perpendicular  to  the  earth-sun  line,  and 
the  vehicle  had  a  spin  period  of  59  seconds.  The  satellite 
encountered  40  day  eclipse  periods  in  both  the  fall  and  the 
spring.  During  these  times,  the  spacecraft  was  shadowed  for 
approximately  1  hour  per  day.  Figure  2  gives  a  schematic 
of  the  orbit.  The  view  is  from  above  the  earth's  north 
pole.  { Re  f . 5 } 

The  environment  of  the  SCATHA  spacecraft  consisted 
of  two  main  regimes,  the  plasmasphere  and  the  plasmasheet. 
The  nature  of  these  regimes,  particularly  the  ion 
composition,  is  important  in  determining  the  magnitude  of 
the  vehicle's  interaction  with  the  environment,  for  example, 
the  sputtering  yield.  The  boundary  between  these  regimes 
(the  plasmapause)  is  well  defined,  but  varies  with  solar 


temperature  is  less  than  1  eV.  The  dominant  ion  is  with 
0*  and  He"^  contributing  approximately  22%  of  the  energy 
density.  The  source  of  these  ions  is  currently  considered 
to  be  the  ionosphere.  {Refs.  6,7,8} 

Overlapping  the  plasmasphere  is  the  ring  current. 
It  consists  of  a  hot,  thin  plasma,  with  densities  normally 
less  than  10®  m“^  and  temperatures  greater  than  10  keV. 
Again,  is  the  dominant  species  (93%  of  the  energy 
density)  while  0"^  and  He"^  are  typically  minor  constituents. 
During  quiet  times  the  principal  source  of  these  * ons 
appears  to  be  the  solar  wind.  During  active  times,  a  large 
0"'  constituent  indicates  ionospheric  sources.  {Refs. 
6,7,3} 

The  plasmasheet  is  the  region  outside  the 
plasmapause,  and  like  the  ring  current,  consists  of  a  hot, 
diffuse  plasma.  Densities  on  the  order  of  10*  -  10®  m“^  and 
temperatures  near  10  keV  .are  common.  At  quiet  times,  the 
composition  is  90%  }{■*■  and  10%  0"^  and  He"^.  At  active  times, 
large  0"^  ratios  (>50%)  are  observed.  This  is  a  mixing  zone, 
containing  particles  of  both  terrestrial  and  solar  origin. 
Large  negative  potentials  are  found  on  spacecraft  in  the 
plasmasheet,  particularly  when  the  satellite  is  in  eclipse. 
Table  1  shows  the  various  ions  in  the  magnetosphere,  and 
their  sources . (Ref s  .  6,7,8} 

An  interesting  feature  of  these  areas  is  the 
variability  of  the  ion  composition  of  the  plasmas  observed 


during  magnetically  active  periods.  In  the  ring  current^  for 
example,  it  has  been  observed  that  0"^  contributes  a  minimum 
of  23%  of  the  storm  time  energy  density,  and  may  contribute 
much  more.  Also,  the  He"^  contribution  increases,  and 
adiabatic  acceleration  of  these  solar  wind  particles  may 
leave  them  with  energies  up  to  32  keV.  Table  2  gives  some 
typical  values  for  storm  time  ring  current  composition. 
{Refs.  7,8} 

The  plasmasheet  region  has  a  direct  source  of  0*,  in 
field  aligned  beams  streaming  from  the  ionosphere.  The 
H‘^/0'^  ratio  of  the  beams  has  been  found  to  be  6/1  {Ref. 

3}.  Additionally,  it  has  been  observed  that  in  the  L  =  6  - 
8  region  (e.g.,  the  SCATHA  orbit),  during  periods  of  high 
magnetic  activity  (Kb.ss5),  the  O'^  density  becomes  comparable 
to  that  of  }{■*•.  (The  Mcllwain  'L'  parameter  is  defined  by 

R=L  cos*0,  where  R  is  the  distance  in  earth  radii,  and  B  is 
the  magnetic  latitude.)  Ks,  is  a  measure  of  the  general 
level  of  magnetic  activity  caused  by  the  solar  wind.  The  Kp 
.scale  is  quas  1 -logar ithmic,  and  runs  from  0-9,  with  0  the 
low  end  and  9  the  high.  On  several  occasions,  the 
density  has  been  found  to  reach  70%.  Table  3  lists  .some 
values  for  the  storm  time  magnetospher ic  bulk  plasma 
composition.  {Refs.  6,7,0} 

3 .  Instruments 


The  data  for  this  experiment  came  primarily  from  the 
University  of  California,  San  Diego  auroral  particles 


experiment,  SC-9.  This  experiment  was  composed  of  five 
detectors;  two  pairs  of  rotating  ion  and  electron  detectors, 
and  one  fixed  ion  detector.  Figure  4  depicts  the  detector 
arrangement.  The  rotating  detectors  scanned  in  orthogonal 
planes  and  were  designated  the  Hi  set  and  Lo  set.  The  Lo 
detector  rotated  from  -20  deg  to  200  deg,  with  70  deg  along 
the  spacecraft  spin  axis.  The  Hi  detector  scanned  the  same 
range,  but  it  was  more  symmetrical,  with  its  midpoint  on  a 
line  parallel  to  the  spin  axis.  Thus,  it  could  not  depress 
more  than  20  degrees  below  the  forward  end  plane.  The  fixed 
detector  looked  radially  away  from  the  spacecraft,  in  the 
same  plane  as  the  Lo  detector.  The  detectors  had  an  angular 
resolution  of  5  by  7  deg.  The  geometry  is  illustrated  in 
Figure  5.  {Ref.  4} 

The  Hi  set  covered  an  energy  spectrum  from  1  eV  to 
81  keV.  The  Lo  rotating  set,  and  the  fixed  detector, 
scanned  from  1  eV  to  2  keV.  For  all  detectors,  a  complete 
energy  scan  required  16  seconds  and  was  covered  in  a  series 
of  64  logarithmic  steps.  The  energy  resolution  at  each  step 
was  approximately  20%.  Additionally,  the  detectors  could  be 
ordered  to  dwell  at  fixed  positions  and/or  energies . {Ref . 4 } 

The  detectors  were  composed  of  three  parts;  an 
electrostatic  analyzer,  an  electrostatic  focussing  lens,  and 
a  spiraltron  particle  sensor.  The  analyzer  provided  the 
energy  differentiation  for  the  system  through  curved  plates 
and  an  applied  voltage.  The  lens  focussed  the  particles  on 


the  sensor  center  by  means  of  two  grids,  one  at  ground  and 
the  *-her  at  the  analyzer  potential.  The  sensor  was  rated 
to  10''  counts  per  second.  {Ref.  4} 

4  .  Ion  Gun 

An  ion  gun  was  installed  on  the  SCATHA  spacecraft  to 
investigate  the  efficiency  of  an  ion  emission  system  in 
modifying  satellite  potentials.  In  particular,  it  could  be 
used  to  develop  negative  voltages  on  the  vehicle.  The 
experiment  utilized  xenon  gas  ionized  by  cathode  discharge 
and  accelerated  by  either  a  1  or  2  kV  potential  drop.  The 
beam  current  could  be  varied  incrementally  from  0.3  mA  to 
2.0  mA.  The  package  also  included  an  electron  source  that 
could  be  configured  to  neutralize  the  particle  beam  or  to 
provide  an  electron  beam.  {Ref.  9} 

C.  OTHER  SATELLITES 

Although  SCATHA  was  the  primary  source  of  data  for  this 
thesis,  two  other  satellites  provided  significant 
information.  These  satellites,  ATS-6  and  ISEE-1,  will  be 
briefly  discussed  in  the  following  paragraphs. 

1.  ATS-6 

ATS-6  was  a  large  three-axis  stabilized  spacecraft 
placed  in  geosynchronous  orbit  in  1974.  Its  major  features 
are  shown  in  Figure  6  .  The  large  antenna  was  significant, 
since  it  fostered  differential  charging  by  shadowing  other 
spacecraft  surfaces.  The  satellite  had  an  environmental 


measurements  package  above  the  spacecraft  dish. 
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Figure  6.  ATS-6  Satellite  {Ref.  12} 


package  contained  rotating  particle  detector  sets  similar  to 
those  found  on  SCATHA.  {Ref. 2} 

2 .  I  SEE 

The  ISEE-1  satellite  was  launched  in  1977  as  part  of 
a  three  spacecraft  mission  to  investigate  the  magnetospher Lc 
plasma.  With  this  in  mind,  the  satellite  was  designed  to 
minimize  absolute  and  differential  charging  by  using 
conducting  materials  and  low  impedance  connections  to 
ground.  Regardless,  the  satellite  was  observed  to  charge  on 
several  occasions.  The  spacecraft  was  equipped  to  measure 
electric  fields,  plasma  waves,  and  plasma  composition.  The 
most  important  of  these,  for  our  purposes,  is  the  plasma 
composition  experiment  which  contained  an  ion  mass 
spectrometer.  (Ret.  10} 

D.  SPACECRAFT  CHARGING 

The  theory  of  spacecraft  charging  will  now  be  reviewed, 
as  differential  charging  will  play  a  significant  role  in 
later  discussions.  In  general,  a  probe  immersed  in  a  plasma 
will  develop  some  potential  relative  to  that  medium.  This 
"floating"  potential  is  determined  by  the  balance  of  the 
currents  incident  upon  the  probe.  The  primary  currents  in 
deciding  the  potential  are;  the  ambient  plasma  electron  and 
ion  flux,  the  photoelectron  emission  due  to  sunlight,  and 
the  fluxes  due  to  secondary  and  backscattered  electrons. 
These  currents  are  schem-at  ica  1  ly  illustrated  in  Figure  7. 
The  left  half  of  the  figure  depicts  typical  eclipse 


conditions  (negative  satellite),  and  the  right  half  the 
currents  present  in  sunlight.  {Refs.  11,12,13} 

In  daylight,  the  dominant  current  is  that  of  the 
photoelectrons  which,  although  material  dependent,  has 
typical  saturation  current  densities  in  the  range  5-50 
}iA/m2  .  This  current  density  overshadows  that  of  the  ambient 
plasma,  which  is  generally  below  1  ^iA/rn^  for  electrons,  and 
.1  JiA/m^  for  ions.  The  value  of  the  secondary  electron 
emission  current  is  highly  dependent  on  both  the  material 


and  the  electron  temperature. 


It  is  thi; 


material 


dependence  which  allows  differential  charging,  even  in 
eclipse.  A  normal  secondary  emission  curve  will  have  a 
yield  greater  than  one  in  some  middle  energy  range 
(typically  50-500  eV),  and  yields  less  than  one  at  higher 
and  lower  energies.  The  point  in  the  high  energy  range 
where  the  yield  decreases  to  one  is  called  the  cross-over 
point  and  is  critical  in  determining  the  effect  of  the 
secondary  emission  current.  The  energy  at  which  this  point 
occurs,  and  the  peak  yield,  varies  with  target  composition. 
Low  energy  ambient  distributions  (T<5  keV}  impacting  the 
spacecraft  will  create  large  secondary  fluxes,  in  fact 
greater  than  the  ambient  electron  flux,  and  the  spacecraft 
will  charge  positively.  High  energy  electron  distributions 
(T^S  keV)  will  produce  lower  secondary  emission  yields,  and 
the  vehicle  will  gain  negative  charge.  Thus,  the  potential 
of  the  satellite  is  determined  by  the  high  energy  electron 


flux.  The  threshold  electron  temperature  for  which  charging 
may  occur  is  typically  in  the  5-10  keV  region.  There  is  an 
additional  citeria,  the  critical  energy,  which  also  must  be 
met.  This  is  an  upper  energy  bound  on  the  distribution 
function  resulting  from  magnetospher ic  convection  processes. 
It  is  also  termed  the  Alfven  boundary,  and  must  exceed  15-20 
keV  for  charging  to  occur.  {Refs.  11,12,13,14,15,16} 

A  simplistic  equation  to  represent  the  situation  may  be 
given  by; 

(1  —  1)  —  lAcitb  ^  “I  m  ma  Ipl-io^o 

Thus,  during  daylight,  when  the  photoelectron  current 
dominates,  the  net  current  is  positive  (toward  the 
spacecraft),  and  the  spacecraft  will  start  to  charge 
positively.  As  this  happens,  the  flow  of  ions  to,  and 
electrons  from,  the  spacecraft  will  be  inhibited,  and  the 
value  of  In.tt  will  decrease.  When  the  value  reaches  zero, 
the  currents  are  balanced  and  the  spacecraft  will  remain  at 
that  potential  until  the  balance  is  disturbed.  For  eclipse 
events,  Iginoeo  is  zero,  and  the  balance  is  effectively 
between  the  ambient  and  secondary  electrons.  For  a  low 
energy  electron  distribution,  the  potential  is  restricted  to 
low  values  (+2  to  -1-5  V),  while  for  high  energies,  t);.' 
poto.ntial  may  reach  -1  to  -20  kV.{Rcfs.  2,11} 

In  addition  to  absolute  charging,  a  spacecraft  may 
charge  differentially.  That  is,  separate  insulated  portions 
of  the  spacecraft  surface  may  charge  to  different 


potentials.  If  two  insulators  are  exposed  to  the  same 


environment,  but  have  different  secondary  emission 
characteristics,  then  their  ability  to  shed  the  excess 
charge  in  the  form  of  a  current  will  differ,  and  a  potential 
will  develop  between  them.  Also,  if  two  portions  are 
exposed  to  differing  environments,  they  will  each  charge  to 
different  potentials.  The  environment  may  differ  over  the 
range  of  a  spacecraft  in  a  number  of  ways.  Frequently, 
pitch  angle  anisotropies  in  the  ambient  plasma  may  cause 
different  fluxes  to  impact  separate  spacecraft  regions. 
Additionally,  if  one  portion  of  the  spacecraft  shadows 
another  continually  (see  ATS-6),  the  shadowed  area  will  not 
emit  photoelectrons  and  will  charge  negatively  compared  to 
the  rest  of  the  satellite.  Also,  in  lower  orbits,  an 
electric  field  induced  in  the  satellite  frame  of  reference 
{E=vxB)  may  cause  variations  in  the  ambient  particle 
distributions.  {Refs.  2,12,13) 

Differential  charging  is  believed  to  be  the  cause  of 
numerous  control  logic  upsets  on  various  spacecraft.  As  the 
differential  potential  increases,  there  will  be  a  tendency 
for  it  to  discharge  in  the  form  of  an  arc.  This  arc  creates 
an  electromagnetic  pulse  which  is  read  by  the  control 
circuits  as  a  command  signal.  For  large  discharges  there 
may  even  be  physical  damage  or  destruction  to  the  arcing 
areas.  {Refs.  17,18) 


Spacecraft  charging  may  be  controlled  by  active  and 
passive  measures.  Active  trials  using  ion  and  electron  guns 
were  carried  out  on  ATS-6  and  SCATHA  with  some  success. 
Passive  attempts,  such  as  those  on  the  ISEE  satellites, 
using  low  impedance  connections  and  conducting  surfaces  were 
also  moderately  successful.  {Refs.  10,19,20} 

E.  SPUTTERING 

The  general  theory  of  sputtering  is  very  significant  to 
this  thesis,  as  it  is  postulated  to  be  the  source  of  the 
observed  spacecraft  generated  ions.  Sputtering  is  defined 
as  the  ejection  of  material  from  solid  surfaces  under  ion 
bombardment.  The  observed  features  of  the  sputtering  yield 
(number  of  target  atoms  sputtered  per  incident  ion)  are;  a) 
its  dependence  on  atomic  number,  b)  its  correlation  with  the 
periodic  table,  c)  its  correlation  with  the  heat  of 
sublimation  of  solids,  d)  its  relation  to  the  efficiency  of 
momentum  transfer,  e)  and  its  dependence  on  incident  ion 
energy,  such  that  it  increased  to  a  maximum  and  then 
decreased.  Figure  8  is  a  plot  of  the  sputtering  yield  for 
aluminum  as  a  function  of  incident  ion  energy  and  mass  {Ref. 
21}.  We  see  that  the  yield  maximum  increases  with  ion  mass, 
and  shifts  to  higher  energies.  These  observations  led  to 
theories  based  on  atomic  collision  cascades.  Such  cascades 
are  created  by  one  incident  ion  impinging  on  the  lattice  and 
transferring  energy  to  other  atoms,  which  in  turn 
participate  in  collisions.  These  theories  are  based  on  the 


and 


simple  elastic  collisions  of  classical  mechanics. 
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particles  to  reach  the  surface  decreases.  {Refs.  21,22,23} 

The  most  comprehensive,  detailed,  and  successful  of  the 
collision  cascade  theories  is  the  Sigmund  -  Thompson  theory. 
This  theory  is  actually  two  theories  developed  separately, 
but  they  complement  one  another  by  explaining  different 
facets  of  sputtering.  The  portion  developed  by  Sigmund 
treats  the  sputtering  yield  of  an  amorphous  target. 
Thompson's  theory  explains  the  energy  distribution  of  the 
sputtered  atoms.  The  development  of  both  these  theories 
will  be  briefly  discussed,  followed  by  a  comparison  with 
experimental  results.  {Refs.  22,23,24,25} 

1 .  The  Sigmund  Theory 

The  Sigmund  theory  focusses  on  the  amount  of  energy 
deposited  in  the  surface  layers  of  the  solid  as  the  driving 
factor  for  the  yield.  The  basic  expression  for  the 
sputtering  yield  is  then  .S  =  rF,  where  r  is  a  constant 
depending  only  on  the  properties  and  state  of  the  target, 
and  F  is  a  function  describing  the  interaction  of  the  two 
particles.  In  order  to  evaluate  this  expression,  Sigmund 


chose  as  his  starting  point  the  sputtering  of  an  amorphous 
target,  and  used  transport  theory  to  describe  the  collision 


cascade  in  random  media.  {Refs.  23,24} 

The  theory  assumes  a  planar  surface  and  an  atom 
starting  its  motion  at  time  t  =  0  and  position  x  =  0.  The 
initial  quantity  of  interest  is  G  ( x,  t  )d^Va  dx.  This 
is  the  average  number  of  atoms  moving  at  time  t  in  a  layer 
(x,dx)  with  a  velocity  (^0,d^V0)  due  to  an  impact  by  a 
single  atom  at  velocity  v.  This  quantity  is  then  used  in  a 
Boltzmann  transport  equation  that  equates  the  initial  number 
G  with  a  later  G  which  is  dependent  on  the  probability  of 
collision.  The  equation  is  then  integrated  to  eliminate 
some  of  the  variables.  Two  functions  are  defined  to 
simplify  the  integration. 


(1-2)  F{x,^0,v)  =  J  G(x,'V0,^,t)dt 
(1-3)  H(x,v)  =  J  F  (  X,  V0,  V)  I  V0,.  I  d^V0 

The  function  F  is  the  total  number  of  atoms  which  penetrate 
the  plane  x  with  a  velocity  (^0,d^V0)  during  the  collision 
cascade  development.  H(x,v)  is  the  sputtering  yield  at  the 
plane  x  for  a  source  of  sputtered  particles  at  x  =  0.  This 
function  may  be  defined  for  either  forward  or  backward 
sputtering.  For  our  purposes,  we  will  choose  backward 


sputtering.  Backward  sputtering  is  emission  directed 
opposite  to  the  initial  direction  of  inc idence  .  f Re f s  .  23,  24  ) 


The  remaining  expression  in  H  is  expanded  in  terms 
of  Legendre  polynomials  and  then  transformed  into  a  moment 
equation  by  multiplying  each  term  by  x"  and  integrating  over 
X.  With  the  assumption  that  the  electronic  losses  to 
electrons  are  small  and  isolated  from  nuclear  collisions, 
the  electronic  and  collisional  cross  sections  may  be 
separated  {Ref. 26}.  This  yields  a  final  expression  which 
includes  the  electronic  stopping  cross  section,  and  the 
differential  cross  section  for  elastic  scattering.  {Refs. 
23,24} 

The  next  step  is  to  determine  the  functions 
describing  the  electronic  and  elastic  cross  sections.  The 
electronic  cross  section  used  is  S.  =  K  S'*,  where  the 
constant  K  is  determined  using  Thomas-Fermi  arguments 
{Ref. 27}.  The  elastic  cross  section  is  derived  from 
classical  mechanics,  and  a  power  approximation  of  the 
Thomas-Fermi  interatomic  potential.  The  equation  is  then 
given  by; 

(1-4)  dCT(T)  =  C  e-"‘T-'^— "  dT 

where  C  is  a  constant,  T  is  the  transferred  energy,  and  E 
the  initial  energy.  The  variable  m  comes  from  the  power 
approximation.  m=l  corresponds  to  Rutherford  scattering,  and 
m=0  approximates  scattering  from  a  Born-Mayer  potential.  If 
these  values  are  substituted  in  the  expression  from  the 
previous  paragraph,  the  following  is  achieved; 
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(1-5)  H(x,E,3)  =  3_  F(x.S.I3) 

4Tt=  N  Cb  Ub 


where  (3  is  the  direction  cosine  of  the  ejection  vector,  and 
F(x,E,|3)  is  a  function  describing  the  distribution  of  the 
energy  deposited  in  the  solid.  Cb  is  a  constant,  and  Ub  is 
the  binding  energy.  In  general,  the  function  F(x,E,(3)  is 
found  to  obey  the  following  relation; 


(1-6) 


F(x,E,(3)  =  a  N  S„(E) 


where  a  is  a  dimensionless  constant  dependent  on  the  ratio 
M2/M1,  N  is  the  atom  density,  and  Sr>(E)  is  the  nuclear 
stopping  power,  given  by  Sr«(E)  “  J  ‘  Additionally,  we 
now  restrict  ourselves  to  sputtering  in  the  unit  normal 
direction.  The  function  a  is  illustrated  in  Figure 
9  .  {Refs. 23, 24  } 

Using  all  of  the  previous  relations,  we  can  then 
define  some  general  expressions  for  the  sputtering  yield. 
For  low  energies  (E<lkeV),  we  have 


(1-7)  S(E,0)  =  3_  a  4M-.  E 

4ti  2  ( M1+M2  )  2  Ub 


Sn(E)  in  this  equation  was  evaluated  by  integrating  da  with 
a  value  of  m=0.  The  Born-Mayer  potential  is  more  accurate 
at  lower  energies  than  the  Thomas-Fermi  appr ox  imat i on .  The 
expression  for  higher  energies  is; 


(1-8) 


S(E,  (3 )  =  3.56  a 


_ tL _ 

(M1+M2 ) 


The  function  Sn{c)  is  the  reduced  stopping  power.  This  is  a 
universal  function  for  Thomas-Fermi  interactions  {Refs. 
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24,27}.  The  variable  E  is  the  reduced  energy,  defined  by 

E  = _ aM:;E _  ,  where  Zi  and  2=  are  the  atomic 

(ZxZ^e^ (M1+M2) ) 

number,  and  a  is  the  Lindhard  screening  radius.  Tabulated 
values  for  Sn(e)  are  shown  in  Table  4.  The  sputtering  yield 
dependence  on  (3  was  roughly  determined  to  be  (cos(0))~^. 
This  dependence  was  derived  empirically  by  Sigmund. 
(Refs . 22,23,24,27} 


The  Thompson  Theory  was  a  complete  theory  of 
sputtering  published  in  1968.  The  portion  involving  the 
determination  of  the  sputtering  yield  was  overshadowed  by 
the  Sigmund  Theory,  but  the  energy  spectrum  determination 
has  endured.  Thompson  started  his  derivation  by  assuming  an 
infinite  solid  with  recoils  at  energy  £=  from  some  primary 
event.  Then,  with  a  density  of  recoil  atoms  g(E=)dE2  per 
unit  time,  each  generating  their  own  collision  cascade,  the 
total  number  of  atoms  slowing  down  through  some  energy  range 
E '  is 


(1-9  ) 


n(E' }  =  p-  q(E=)  U(E=,E' )dE: 


The  mean  rate  of  energy  loss  is  given  by  vdE'/dx,  where  v  is 
the  velocity.  Assuming  an  isotropic  distribution,  one  may 
define 


( 1-10) 


n(E'  ,f )dE'dft’  =  p- 


p-  Q(E^)U(E: 
Jb-v  dE'/dx 


)  dEs  dE'  dfl' 
4Tt 


From  this,  the  flux  may  be  determined, 

(1-11)  <J)  ( E  '  ,  r  )  dE  '  dii '  =  V  n  (  E  '  ,r  )  co39dE  '  ' 

The  theory  now  assumes  that  the  infinite  solid  is 
cut  in  half  and  the  flux  is  observed.  Also,  some 
expressions  are  needed  for  the  functions  above.  jKEaiE')  is 
the  number  of  displaced  atoms  and  may  be  approximated  by 
[BEa/E',  where  (3  is  a  constant  of  order  unity.  Similarly, 
dE'/dx  may  be  replaced  with  E’/D,  where  D  is  the  interatomic 
spacing.  Using  these  relations  we  obtain 

(1-12)  <i>(E' ,r  )dE'dft’=13_^  j'“q(E=  )E  =  dE2  cos9  dii'  dE' 

This  function  i  is  the  flux  inside  the  surface,  and  must  be 
transformed  to  that  observed  outside.  This  is  accomplished 
by  assuming  a  binding  force  normal  to  the  surface  which  will 
decrease  the  normal  kinetic  energy,  but  leave  the  parallel 
component  unaffected.  Thus,  the  effect  will  be  a  bending  of 
the  particles  trajectory.  With  the  energy  equations,  and 
the  relationship  between  the  inner  and  outer  angles  in  hand, 
we  come  to  the  final  expression 

(1-13)  'f'(E.0)dEd5l=  (3  D  cos0  P”  q(E2)E2dE3  dEdn 

4ti(  H-Ub/E)  Ja¬ 

in  this  equation,  Ub  is  the  binding  energy,  E  is  the  emitted 
particle  energy,  and  *  is  the  angle  of  the  trajectory  from 
the  surface  normal.  The  chief  features  are  the 


proportionality  to  E/(E+Ub)^  and  the  cos0  dependence. 


Figure  10  is  a  plot  of  this  function  versus  energy.  For 
high  energies,  the  spectrum  will  fall  off  as  l/E^,  and  for 
low  energies,  the  binding  energy  will  have  an  impact.  The 
maximum  of  this  curve  occurs  at  U0/2.  Figure  11  is  a  plot 
of  the  measured  energy  d istr ibut ion, with  the  model  overlaid. 
It  can  be  seen  that  up  to  approximately  1000  eV,  the  model 
matches  the  observed  spectrum  quite  closely.  Since  the 
majority  of  the  sputtered  particles  are  in  the  energy  region 
that  follows  the  1/E^  curve,  this  theory  is  quite  effective 
in  describing  the  data.  {Refs.  22,24} 


To  evaluate  the  applicability  of  these  theories,  we 
must  determine  how  closely  they  recreate  experimental 
results.  Therefore,  we  will  look  at  the  sputtering  yield  as 
a  function  of  ion  energy,  angle  of  incidence,  incident  and 
target  masses,  and  surface  binding  energy.  The  energy 
dependence  of  the  sputtering  yield  is  contained  in  the 
energy  dependence  of  the  stopping  power.  Calculated  yields 
using  equations  1-7  and  1-8  are  compared  with  experimental 
values  in  Figure  12  for  copper  (Ref.  23}.  The  solid  line 
is  the  high  energy  result,  and  the  dotted  line  is  the  plot 
for  the  low  energy  equation.  The  experimental  results  are 
marked  by  the  various  geometric  shapes.  In  addition, 
calculated  and  measured  yields  for  various  materials  are 
listed  in  Table  5  (Ref.  23}.  It  can  be  seen  that  the 
correlation  is  reasonable,  normally  within  a  factor  of  two. 


(eV) 


Some  of  the  experimental  yields  are  questionable,  as  the 
sputtering  yield  is  known  to  have  a  dose  dependence.  The 
Sigmund  Theory  also  tends  to  overestimate  sputtering  at  very 
low  energies  (<100  eV).  In  most  cases,  experimental  data 
indicates  that  the  yield  goes  to  zero  in  the  50-100  eV 
region . {Ref s .  22,23,24} 

The  sputtering  yield  is  a  function  of  Mi  directly 
and  in  the  factor  oi(M2/Mi).  Due  to  the  above  mentioned  dose 
effects,  yield  data  for  masses  are  generally  normalized  to 
the  self  sputtering  yield  or  the  argon  sputtering  yield. 
Data  for  various  targets  are  presented  in  Figure  13  (Ref. 
23}.  The  top  plot  displays  results  for  impacts  on  silicon. 
The  mass  of  the  sputtering  ion  in  AMU  is  along  the 
horizontal  axis  and  the  yield  of  the  ions  normalized  to  the 
argon  sputtering  yield  on  the  vertical  axis.  The  other 
plots  are  similar  for  copper  and  gold.  The  agreement  with 
theory  for  silicon  is  quite  good,  but  decreases  for  the 
heavier  targets.  Sputtering  yields  for  and  0"^  on  silicon 
are  illustrated  in  Figure  14.  Since  in  general  we  will  be 
discussing  silicon  or  silicon  dioxide  as  a  target,  we  will 
disregard  the  lack  of  agreement  for  heavier  targets  in  the 
model  which  we  develop.  {Refs.  22,23,24} 

The  angular  dependence  of  the  sputtering  yield  is 
considered  to  be  proportional  to  (l/cosS)".  The  exponent, 
n,  is  a  variable  quantity  ranging  from  a  little  less  than  1 
to  1.7.  For  the  situations  we  will  be  discussing. 


with 


M=/Mi  >  1,  the  value  of  n  is  generally  taken  to  be  1.7. 
Figure  15  is  a  plot  of  the  dependence  on  angular  incidence 
from  experimental  data  {Ref.  23}.  The  dependence  ranges  up 
to  an  incident  angle  of  70  degrees.  Beyond  this  point,  the 
sputtering  yield  decreases  rapidly  to  0.  If  we  attempt  to 
model  using  only  the  0-70  degree  range,  we  will  probably  be 
underestimating  the  yield  by  approximately  20°s.  {Refs. 
22,23,24} 

The  yield  depends  directly  on  the  binding  energy, 
thus  the  determination  of  this  energy  is  quite  important. 
Sigmund  suggested  that  the  appropriate  value  for  the  binding 
energy  would  be  the  sublimation  energy  of  the  material. 
This  generally  leads  to  sputtering  yield  values  that  are 
slightly  high.  Recent  work  has  given  theoretical  backing 
for  a  higher  value,  but  no  definite  numerical  quantity  was 
put  forward  {Ref.  28}.  It  has  been  suggested  that  the  value 
providing  the  best  results  is  17/10  of  the  sublimation 
energy  {Ref.  29}.  If  this  value  is  used  to  recalculate  the 
numbers  in  Table  5,  we  get  the  results  in  Table  6.  It  is 
obvious  that  in  general,  the  calculated  yields  in  Table  6 
are  much  closer  to  the  experimental  yields,  than  those 
calculated  for  Table  5.  For  calculations  in  this  thesis,  we 
will  use  the  larger  value  of  the  binding  energy.  {Refs. 
22,23,24} 
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Angle  of  Incidence 


Angle  of  Incidence 


Figure  L5.  The  Dependence  of  the  Sputtering  Yield 
on  Angle  of  Incidence  {Ref. 23} 
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The  theory  was  devised  for  amorphous  monatomic 


targets  bombarded  by  atomic  ions.  It  is  not  so  readily 


applicable  to  the  bombardment  of  molecular  solids,  as  will 


be  necessary  in  later  calculations.  Since  several  of  the 


reference  data  points  are  concerned  with  bombardment  vich 


molecular  oxygen,  it  will  be  relevant  to  discuss  this  type 


of  impact.  The  case  of  molecular  ions  on  molecular  targets 


may  be  treated  through  computer  simulation,  but  it  is  an 


extensive,  complicated  project,  as  the  probabilities  for 


each  collision  combination  must  be  calculated.  In  most 


cases,  it  is  more  appropriate  to  utilize  experimental  data 


if  it  is  available.  Such  data  is  often  difficult  to 


obtain.  Only  recently  has  an  effort  begun  to  systematically 


tabulate  the  known  experimental  results  and  compare  them 


with  the  theory.  The  presently  tabulated  data  concerns 


atomic  ions  on  monatomic  targets  {Ref.  30}.  A  great  deal  of 


data  is  also  distributed  randomly  through  the  literature  on 


other  ion-target  combinations. 


An  additional  bit  of  information  is  available  for 


the  yields  of  monatomic  molecules.  It  has  been  found  for 


several  atoms  that  the  sputtering  yield  per  atom  for 


molecular  ions  is  slightly  higher  than  that  for  atomic  ion: 


(Ref.  31}.  The  results  are  listed  in  Table  7.  The  yield: 


are  greater  by  approximately  a  factor  of  two,  while  for 


silicon  they  range  from  1.15  to  1.30.  There  i: 


difference,  but  it  is  not  drastic.  While  it  is  difficult  to 
accurately  model  the  more  complicated  situations,  it  should 
be  possible  to  obtain  some  approximate  results  for  monatomic 
molecules  that  are  within  a  factor  of  two  of  the  actual 
results . 


F.  NASA  CHARGING  ANALYZER  PROGRAM 

The  NASA  Charging  Analyzer  Program  (NASCAP)  was 
developed  to  accurately  model  the  dynamics  of  spacecraft 
response  to  realistic  plasma  environments.  NASCAP  can 
effectively  simulate  the  charging,  in  both  laboratory  and 
magnetosphere  environments  of  objects  that  are 
geometrically,  materially,  and  electrically  complex.  The 
program  utilizes  a  timestep  procedure,  calculating  quasi¬ 


static  steps 


The  dynamics  are  driven  by  charge 


accumulation  from  external  sources,  charge  depletion,  and 
conduction  in  dielectrics.  Each  timestep  includes  a  full 
three  dimensional  electrostatic  potential  calculation.  The 
calculation  of  Poisson's  equation  over  a  grid  is  time- 
staggered  with  a  procedure  in  which  incident  charged 
particle  fluxes,  leakage  currents,  emission  currents,  and 
induced  spacecharge  effects  are  computed  based  on  the 
current  quasistatic  conditions.  NASCAP  is  limited  to 
situations  where  the  Debye  length  is  long  compared  to  the 
dimensions  of  the  examined  object.  This  condition  is  met 
for  the  SCATHA  spacecraft.  NASCAP  can  output  a  time  history 


of  spacecraft  charging,  potential  contours,  charge  contours, 
cur-rent  contours,  and  particle  trajectories. 


A.  PREVIOUS  OBSERVATIONS 


Data  taken  from  three  satellites,  AT3-5,  ATS-6,  and 
ISEE-i,  Will  be  discussed,  prior  to  presenting  P78-2 
(3CATHA)  results.  Observations  of  the  three  primary 
phenomena  (triangle  peak,  shadow  peak,  diffuse  background; 
on  satellites  prior  to  SCATHA  will  be  illustrated. 

1.  ATS-5 


The  presence  of 

ions  below 

the  charging 

peak  was 

first  reported 

for  ATS- 

-5  {Ref.l}. 

ATS-5  data 

taken  in 

eclipse  in  1969 

and  1970 

showed  ion 

fluxes  in  a  broad  range 

of  energies  below  the  charging  peak.  A  typical  example  of 
such  data  is  shown  in  Figure  16.  The  data  is  presented  in 
grey  scale  spectrogram  format,  a  common  method  of  displaying 
the  data  concisely.  The  figure  is  vertically  separated  into 
three  components.  The  top  primarily  displays  magnetic  field 
information,  unnecessary  for  our  purposes.  The  remaining 
two  sections  indicate  the  count  rates  of  electrons  and  ions 
in  the  detector  versus  time  and  energy.  The  horizontal  axis 
defines  increasing  time,  and  in  this  case  is  labelled  with 
the  hours  of  day  274  of  1970.  The  vertical  axi.s  is 
associated  with  the  energy  of  the  particles.  The  minimum 
(~50  eV)  of  the  energy  axis  is  located  between  the  two 
sections,  and  energy  increases  upward  for  electrons,  and 


downward  for  ions  to  a  maximum  of  50  keV.  The  magnitude  of 
the  detected  countrates  is  indicated  by  the  shading  of  the 
figure  from  black  to  white,  with  black  cor respond i ng  to  zerc 
counts,  and  white  with  maximum.  If  the  detector  is 
saturated  at  a  certain  time  in  an  energy  channel,  the 
spectrogram  will  contain  a  black  spot  at  that  time  and 
energy.  A  feature  such  as  this  can  usually  be  distinguished 
from  a  zero  count  rate  because  saturation  points  are 
normally  surrounded  by  intense  count  regions,  revealed  as 
bright  white  areas. 

The  data  displayed  in  Figure  16  are  from  day  274  of 
1970.  The  satellite  enters  a  region  of  hot  plasma  at  0450, 
as  indicated  by  the  increased  intensity  of  high  energy 
(>4500  eV)  electrons.  The  satellite  charging  peak  may  be 
seen  as  the  intense  white  region  in  the  proton  spectrum 
between  0620  and  0720.  From  the  energy  of  this  region,  we 
can  infer  a  satellite  potential  of  approximately  .75-2  kV. 
The  black  spots  on  the  peaks  are  an  example  of  how  detector 
saturation  is  displayed.  The  broad  spectrum  of  ions  is 
evident  in  the  energy  bands  below  the  charging  peak  as  the 
grey  area  that  clearly  extends  from  1  kV  to  about  300  V,  and 
may  even  continue  to  lesser  energies.  These  particles  did 
not  appear  during  every  charging  event,  so  the  source  was 
difficult  to  determine  and  no  complete  explanation  was  put 


or ward . 


contamination  of  the  ion  data.  Secondary  or  shadow  peaks 
were  observed  at  energies  below  the  inferred  satellite 
potential.  An  example  of  such  data  is  illustrated  in  Figure 
17.  This  figure  is  a  grey  scale  spectrogram  similar  to  the 
one  shown  in  Figure  16.  The  data  presented  are  from  the 
North-South  detector  on  ATS-6  on  day  243  of  1974. 

This  plot  displays  two  phenomena.  The  first  is  the 
ion  peaks,  which  are  clearly  presented  elsewhere  {Ref.  2}. 
These  ions  were  denoted  "spacecraft  generated  ions".  The 
second  is  an  apparently  analagous  phenomenon  in  the  electron 
spectrum.  These  observations  are  the  so-called  "Minnesota 
spots".  In  addition  to  the  electron  and  ion  spectra,  this 
figure  displays  the  detector  pitch  angle  in  the  top  section, 
which  is  related  to  the  rotation  angle. 

The  spots  are  seen  in  the  electron  spectrum  between 
100  and  300  eV  as  intense  white  points  that  appear  to  form 
triangle  shapes.  If  one  compares  the  periodicity  of  the 
electron  triangles  with  the  detector  look  angle  periodicity, 
a  direct  correspondence  is  evident.  As  the  top  of  the 
spacecraft,  where  the  detectors  are  located,  was  in  shadow 
at  this  time,  the  source  electrons  were  determined  to  be 
secondaries  from  the  surface  of  the  vehicle.  These 
electrons  had  been  emitted  from  a  surface  differentially 
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Figure  17.  ATS-6,  Spectrogram,  Day  248  of  1974  {Ref. 2} 
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rharged  with  respect  to  the  detectors,  and  thus  the 


distribution  had  been  accelerated  by  a  potential  drop  0. 


The  source  of  the  flux  was  eventually  determined 


secondary  emissions  from  the  University  of  Minnesota 


detectors  also  located  on  the  satellite  measurements  box. 


{Ref.  2} 


3 .  I  SEE 


Negative  charging  events  on  the  I.SEE  satellites  were 


■latively  rare,  as  the  satellites  were  specifically 


designed  to  avoid  the  problem  of  spacecraft  charging. 


1  did,  however,  charge  to  significant  negative  levels  on  a 


few  occasions.  During  one  such  event,  analysis  of  the  ion 


mass  composition  revealed  some  interesting  features.  Figure 


13  is  a  presentation  of  the  ISEE-1  mass  spectrometer  reading 


from  March  17,  1973.  The  vertical  axis  is  a  count  rate,  and 


the  horizontal  axis  is  the  applicable  mass  channel.  The 


mass  channel  is  related  to  the  atomic  mass  by. 


(  2-1) 


AMU  =  _ 3515 _ 

Mass  Channe 1 - • ^ 


Thus,  mass  channel  54  corresponds  to  H"^,  15  to  N"^,  and  14  to 


0".  {Ref.  10} 


At  0633  UT  (Universal  Time),  the 


ms  t  r  ume  nt 


responded  with  a  broad  peak  in  the  high  Atomic  Mass  Unit 


'iMU=30-100)  channels.  This  peak  had  a  maximum  in  channel 


10  corresponding  to  32  AMU.  There  was  also  a  peak  a' 


channel  54,  H*.  The  normal  value  of  the  ambient  atomic 


S  61.2’O/SiNnOO 


Figure  18,  ISEE-1,  Mass  Spectrometer,  March  17,  1978 


MASS  CHANNEL' 
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oxygen  is  indicated  by  the  solid  black  line  under  the  peak. 
The  data  from  0721  UT  are  quite  similar,  except  that  at  this 
point  the  maximum  has  shifted  to  channel  11,  28  AMU.  Again, 


this  is  larger  than  the  indicated  ambient  oxygen.  At  the 
time  these  observations  were  first  made  it  was  concluded 
that  these  molecules  were  probably  hydrocarbons  and  other 
outgassing  products.  In  light  of  spacecraft  surface 
sputtering,  however,  one  may  postulate  a  different  flux 
source.  That  is,  these  particles  are  ions  emitted  from  the 
spacecraft  surface.  Since  the  surface  is  primarily  silicon 
dioxide  (SiOs),  the  ejected  particles  would  largely  be  Si, 
0,  and  O2.  The  masses  of  these  are  28,  16,  and  32  AMU 


respectively,  numbers  which  would  fit  nicely  into 
observed  data.  {Ref. 10} 


the 


B.  P78-2  (SCATHA) 

Observations  of  spacecraft  generated  ions  are  presented 
for  five  days.  The  first  two  are  for  eclipse  charging 
events  in  1979,  the  second  pair  are  for  ion  gun  induced 
charging,  and  the  last  is  for  eclipse  charging  in  1981. 

1.  Day  83  of  1979 

This  first  eclipse  event  illustrates  the  triangle 
peak  and  its  associated  angular  dependence.  Figure  10  is 
the  grey  scale  spectrogram  for  a  one  hour  period  of  day  "1 
of  1979.  This  is  a  composite  of  the  electrons  from  the  Hi 
detector,  and  the  ions  from  all  three  detectors.  Each  of 
the  ion  displays  has  zero  energy  at  the  top  and  increases 

4  8 


k  jiCS 


energy  downward.  The  electron  presentation  is  reversed, 
with  low  energies  at  the  bottom,  increasing  upward.  Between 
1740  and  1820,  the  triangle  structures  can  be  seen  in  the  Lo 
detector  data.  From  the  Hi  detector  data,  it  can  be  seen 
that  the  spacecraft  charges  negatively  to  about  200  v  at 
1750,  decreases  to  100  V  at  1815,  and  Increases  rapidly  to 
300  V  at  1825.  Multiple  shadow  peaks  are  visible  during 
this  period  in  both  the  Hi  and  Fix  detector  data. 

Figure  20  is  a  plot  of  the  log  ion  distribution 
function  {phase  space  density)  versus  energy  for  the  Hi 
detector.  The  spacecraft  is  charged  to  -250  V,  in  eclipse, 
as  indicated  by  the  peak  in  phase  space  density  (f)  at  250 
eV.  The  shadow  peak  is  the  secondary  peak  at  78  eV.  At 
this  time,  the  detector  is  parked  at  92  degrees, 
approximately  parallel  to  the  spacecraft  spin  axis. 

Figure  21  is  a  diagram  of  the  spacecraft  potential, 
triangle  peak  energy,  and  Lo  detector  look  angle  versus 
time.  The  energy  of  the  peak  varies  directly  with  the  look 
angle,  with  the  energy  minimum  occurring  at  the  angle 
maximum,  that  is,  looking  down  and  away  from  the  spacecraft 
(where  it  approaches  closest  to  the  spacecraft  surface). 
The  shadow  peak  energy  maintains  a  relatively  constant  ratio 
with  satellite  potential.  In  this  case,  Ep/Et  is  between  6 
and  7,  where  Ep  is  the  satellite  potential,  and  Et  is  the 
triangle  peak  minimum  energy.  This  is  significant,  as  we 


would  expect  a  differentially  charged  surface  to  maintain  a 
constant  ratio  over  time  with  the  satellite  mainframe. 

2.  Day  86  of  1979 

The  second  eclipse  charging  event,  day  36,  is 
similar  to  that  of  day  83.  Figure  22  is  a  grey  scale 
spectrogram  of  the  ion  data  for  a  one  hour  period  on  day  36 
of  1979.  The  triangle  structures  are  again  visible  from 
1545  to  1740  in  the  Lo  detector  display.  From  the  Hi  data 
we  see  that  charging  begins  at  approximately  1640,  and  the 
potential  varies  widely  throughout  the  period.  A 
fluctuation  in  the  satellite  potential  due  to  a  change  in 
the  electron  Alfven  boundary  is  seen  at  1705.  The 
spacecraft  achieves  a  maximum  potential  of  -600  V  at  1725. 
During  the  period,  a  shadow  peak  is  clearly  visible  at 
energies  lower  than  the  charging  peak  in  both  the  Hi  and  Fix 
detector  data. 

Figure  23  plots  the  ion  distribution  from  the  Hi 
detector.  The  satellite  potential  is  -340  V  as  indicated  by 
the  charging  peak.  There  is  a  well  defined  shadow  peak  at 
100  eV.  In  addition,  there  are  various  other  small  peaks 
below  the  charging  peak  that  may  indicate  other  shadow 
peaks.  Frequently,  other  sh.adow  peaks  appear  in  the 
spectrograms,  maintaining  a  constant  separation  ratio 
between  each  other,  and  varying  as  the  main  peak  varies. 


The  detector  is  again  parked  parallel  to  the  spacecraft  spin 
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The  strong  dependence  between  the  triangle  peak 


minimum  energy  and  spacecraft  potential  is  illustrated  by 


Figure  24  for  this  more  variable  day.  It  is  clear  that  the 


minimum  energy  of  the  triangle  peak  increases  with  the 


increasingly  negative  vehicle  potential.  The  ratio,  Ep/Et, 


■s  maintained  between  6  and  6.5.  Again,  the  energy  of  the 


shadow  peak  varies  directly  with  the  detector  angle. 


achieving  minimum  energy  at  the  detectors  maximum  angular 


extent.  The  shadow  peaks  usually  lose  their  definition 


below  110  degrees. 


These  two  examples  illustrate  typical  observations 


during  eclipse  passages  for  which  negative  charging 


occurred.  By  comparison  with  ATS~6,  we  would  expect 


similar  observations  during  daylight  negative  charging 


events.  ATS-G,  however,  was  stabilized  and  had  its  large 


solar  array  to  provide  a  shadowed  area  which  would  simulate 


an  eclipse.  SCATHA  had  no  such  array,  and  in  addition  was 


spinning  so  that  no  area  was  shadowed  long  enough  to  charge 


to  a  high  level.  Observations  of  daylight  charging  showed 


that  the  spacecraft  potential  was  highly  spin  modulated. 


making  it  difficult  to  use  these  data  for  this  project 


However,  induced  charging  events  with  the  ion  gun  provide  an 


additional  range  of  daylight  data.  Two  examples  of  such 


experiments  are  presented. 


msa 


PEAK  ENERGY  (eV) 


The  first  example  of  ion  gun  induced  charging  is  day 


94  of  1979 

The 

ion  gun  was  in  operation  at 

1  kV  and  1 

mA  . 

Figure  25 

is  the 

grey  scale  spectrogram 

for  a 

1  hour  period 

from 

that 

day . 

The  data  presented 

are 

from  two 

X  l  i 

detectors. 

w  i  th 

the  Lo  detector  in  the 

upper 

d iagram. 

and 

the 

Fix 

in  the 

lower.  The  ion  gun 

1  3 

turned  on 

at 

approximately  1405  and  the  spacecraft  charges  promptly  to  a 
final  value  of  -60  V.  Shadow  peaks  occur  at  11  and  20  eV  in 
the  Lo  detector  between  1420  and  1430,  and  between  1440  and 
1455  . 

A  curious  feature  is  the  appearance  of  twin  shadow 
peaks  in  the  Lo  data  at  1440,  with  bands  between  in  which  no 
ions  were  detected.  When  this  phenomena  began,  the  detector 
had  just  been  parked  at  an  angle  of  20  degrees  to  the 
spacecraft  spin  axis,  away  from  the  satellite  body.  The 
detector  had  previously  been  scanning.  When  the  data 
printouts  are  examined,  three  bands  clearly  emerge;  one  at 
11  eV,  another  at  22  eV,  and  the  last  at  30  eV.  We  may 
tiiat  certain  energies,  or  trajectories,  do  not  reach 
the  detector  at  that  angle.  These  events  indicate  something 
about  the  location  of  the  emitting  surfaces.  From  the 
striations  in  the  Lo  data,  we  would  expect  non-emitting  or 
non -d  i  f  f  er  en  t  ia  1  ly  charged  surfaces  to  be  sandwiched  betwecnn 
emitting  regions,  thus  giving  the  dark  regions  where  no  ions 
■sre  detected.  Importantly,  the  Ki  detector  is  viewing  in  a 


ASHaNa 


(A»)  ADaana 


Figure  25.  SCATHA,  Spectrogram,  Day  94  of  1979 


similar  direction  at  this  time,  and  gives  a  similar 
response.  The  Fix  detector,  viewing  70  degrees  from  the  Lo 
and  Hi,  reports  a  completely  different  environment.  It 
does,  however,  show  a  similarity  in  that  a  large  energy  gap 
exists.  This  indicates  mixed  differentially  and  non- 
differentially  charged  surfaces  as  above. 

4 .  Day  200  of  1979 

Our  fourth  example  is  during  an  ion  gun  operation  at 
1  kV,  1  mA .  Figure  26  is  a  grey  scale  spectrogram  for  a  two 
hour  timeframe.  The  data  are  taken  from  the  Hi  detector. 
In  this  figure,  the  spacecraft  is  charged  to  approximately  - 
700  V,  except  during  the  two  five  minute  periods  when  the 
ion  gun  was  in  "trickle”  mode.  Trickle  mode  is  defined  as 
the  gun  in  operation  with  the  accelerating  voltage  off,  but 
current  still  applied  to  the  diode.  Thus,  the  generated 
ions  leak  out  at  low  energies.  Triangle  peaks  are  clearly 
visible  in  the  ion  display.  This  day  is  one  of  the  few  in 
which  these  peaks  become  evident  in  the  Hi  detector.  The 
ratio,  Ep/Et,  for  these  figures  is  in  the  9-10  range, 
indicating  that  the  Hi  detector  is  looking  at  an  emission 
point  different  from  that  seen  by  the  Lo  detector  on  days  83 
and  86  (or  different  processes  are  at  work). 

Figure  27  shows  the  distribution  functions  when  the 
Lo  and  Fixed  detectors  are  looking  approximately  parallel, 
radially  away  from  the  spacecraft.  The  Hi  detector  is  near 
its  maximum  rotational  angle.  The  satellite  potential  is 
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approximately  -700  V.  The  Hi  and  Fix  detectors  exhibit  a 
broad  array  of  low  energy  ions,  but  such  a  spectrum  is  not 
evident  in  the  Lo  detector.  The  Lo  detector  does  show  a 
sharply  defined  shadow  peak.  The  data  have  a  similar  form 
throughout  this  period,  with  subtle  variations  in  the 
amplitude  and  energy  of  the  diffuse  spectrum,  particularly 
with  angle. 

5.  Day  92  of  1981 

Figure  28  is  the  grey  scale  spectrogram  from  a  1981 
eclipse  charging  event.  The  figure  displays  4  hours  of  data 
from  the  Hi  detector  for  both  ions  and  electrons.  The 
eclipse  charging  event  is  visible  in  the  illustration 
between  0810  and  0945  UT,  with  potentials  varying  from  -100 
to  -400  V.  During  this  period,  a  definite  shadow  peak  is 
also  evident,  in  the  -60  to  -150  V  region. 

Figure  29  is  a  plot  of  the  ion  distribution 
functions  for  day  92.  The  log  distribution  function  has  a 
charging  peak  at  0820  at  90  eV  with  a  value  of  6.07.  This 
value  is  comparable  to  those  on  days  83  and  86  of  1979. 


Thus  the  detected  countrates  are  as  high  as  those 
encountered  early  in  the  mission.  This  indicates  that  the 
source  of  the  detected  particles  is  not  dependent  on 
satellite  lifetime,  and  hence  that  outgassing  is  not  a  major 
contributor  to  the  observed  satellite  generated  ion 
phenomena . 
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Figure  29.  Day  92  of  1981,  Ion  Distribution  Function 


A.  POSSIBLE  ION  SOURCES 


The  ion  source  must  be  able  to  explain  three  different 
phenomena,  the  triangle  peaks,  the  shadow  or  mirror  peak, 
and  the  diffuse  spectrum.  The  prime  feature  that  the 
observations  have  in  common  is  the  appearance  of  the  subject 
ions  at  energies  below  that  provided  by  the  spacecraft 
potential.  Thus,  they  must  be  generated  in  regions  near  the 
spacecraft  where  they  do  not  experience  the  full  plasma- 
satellite  potential.  There  are  two  prime  candidates  for  the 
source,  outgassing  and  sputtering.  Outgassing  is  the 
emission  from  the  satellite  of  contaminant  molecules  and 
atoms  with  energies  near  the  temperature  of  the  satellite. 
These  are  primarily  trapped  hydrocarbons  and  other 
atmospheric  molecules.  The  outgassing  rate  may  be  increased 
by  the  inclusion  of  efflux  from  the  attitude  control  jets, 
and  the  ion  beam  system.  Sputtering  is  the  emission  of 
charged  and  neutral  atoms  from  the  satellite  surface  due  to 
ion  bombardment,  and  its  theory  has  been  described 
previously.  A  qualitative  discussion  of  both  possibilities 
follows,  followed  by  a  numerical  calculation  of  the  number 
of  ions  that  would  reach  the  detectors  due  to  Sigmund- 
Thompson  sputtering  from  differentially  charged  surfaces. 


As  the  contaminant  particles  leave  the  spacecraft 
they  are  subject  to  ionization  from  photons,  and  collisions 
with  ambient  ions  and  electrons.  As  the  molecules  are 
ionized  at  various  distances  from  the  spacecraft  they  will 
return  in  a  broad  range  of  energies.  This  is  a  possible 
explanation  of  the  observed  diffuse  spectrum.  The  question 
is  then,  how  many  particles  outgas,  and  what  percentage  of 
these  are  ionized?  The  problem  of  outgassing  and  the 
outgassing  yield  has  been  studied  previously  {Refs.  32,33}. 
The  dominant  mechanism  for  the  ionization  of  the  outgassing 
molecules  has  been  determined  in  sunlight  to  be 

photoionization,  and  in  eclipse  to  be  collision  with 

electrons.  The  electron  ionization  rate  is  approximately 
two  orders  of  magnitude  less  than  the  photoionization  rate. 
In  eclipse,  however,  the  photoionization  rate  goes  to  zero 
and  the  electron  collisions  become  dominant.  In  all 
probability,  the  electron  ionization  rate  will  increase 
slightly  in  eclipse,  as  the  electron  density  will  increase 
due  to  space  charge  effects.  The  ionization  rate,  however, 
will  not  exceed  that  of  the  photoionization.  (Refs.  2,32) 

The  distance  within  which  the  ions  may  be  returned 
to  the  spacecraft  will  be  the  width  of  the  surrounding 
plasma  sheath.  The  sheath  distance  at  geosynchronous  orbit 
is  approximately  10Pt  m.  With  this  information  and  the 
ionization  rates  above,  it  may  be  shown  that,  at  a  maximum. 


approximately  .02%  of  outgassing  silicon  atoms  will  be 


ionized.  In  a  central  force  field,  with  conservation  of 
angular  momentum  and  energy,  calculations  reveal  that  all 
the  particles  ionized  in  the  sheath  will  be  returned  to  the 
spacecraft.  The  final  result  then,  is  that  .02%  of  the 
outgassing  molecules  will  be  reattracted  to  the  satellite. 
{Ref.  32} 

When  the  spacecraft  is  first  launched,  the  ionized 
outgassing  flux  may  be  sufficient  to  provide  a  significant 
number  of  particles  at  the  detectors.  Outgassing,  however, 
is  an  extremely  time  dependent  phenomenon.  Once  the 
particles  leave  there  is  no  way  to  replenish  them,  so  the 
outgassing  rate  must  decrease  with  time.  Various 
experimental  data  have  placed  the  outgassing  as  proportional 
to  t“",  where  n  varies  from  0.8  to  1.6,  or  as  proportional 
to  exp ( -t/41 . 6 ) ,  with  t  in  days  {Refs.  32,33}.  The  longest 
e-folding  time  then  would  be  41.6  days.  Even  this  would 
lead  to  a  decrease  in  the  outgassing  flux  by  three  orders  of 
magnitude  over  the  period  of  a  year.  As  noted  in  day  92  of 
1981  above,  however,  the  shadow  peaJcs  have  not  diminished 
over  time  in  any  significant  way.  This  leads  to  the 
conclusion  that  the  outgassing  flux  may  be  a  contributor, 
but  is  not  the  prime  source  of  the  ions. 

What  effect  do  the  thrusters  and  ion  gun  have  on 
this  conclusion?  An  experiment  was  performed  aboard  the 
SCATHA  spacecraft  to  measure  the  contamination  (impurities 


a!: 


mixed  with  the  hydrazine  fuel)  on  the  satellite  due  to  the 
thruster  operation.  Two  sensors  were  placed  on  the 
satellite,  one  in  the  bellyband,  and  one  on  the  forward  end. 
The  sensors  were  capable  of  detecting  5  ng/cm^  particle 
deposition  and  a  current  as  low  as  10”^*  A.  The  thrusters 
were  located  on  the  aft  end  of  the  vehicle,  and  were  fired 
once  per  week  during  the  bulk  of  the  flight  time  to  correct 
for  precession.  The  results  of  the  experiment  indicated 
that  no  measurable  flux  of  thruster  contaminants  was 
returned  to  the  spacecraft.  {Ref. 34} 

The  ion  gun  is  also  an  unlikely  source  of 
contaminants  to  the  detectors.  As  noted  above,  it  is 
located  on  the  aft  end  of  the  spacecraft  and  directed  away 
from  it.  During  most  normal  modes  of  operation,  the  energy 
and  angular  momentum  of  the  exiting  particles  would  prohibit 
them  from  reaching  the  detectors.  There  is  evidence  that 
some  ions  are  returned  to  the  detectors,  but  they  are 
observed  at  energies  greater  than  the  satellite  potential 
{Ref.  35}.  There  is  a  possibility  that  ions  emitted  during 
trickle  mode,  at  low  energies,  may  be  reattracted,  and 
contribute  to  the  flux  observed  at  the  detectors.  This  is 
probably  only  a  small  intermittent  component,  as  the 


With  just  the  emissions  of  the  contaminant  flux,  it 
is  still  numerically  conceivable  that  the  outgassing  flux 
could  produce  sufficient  particles  to  account  for  the 
observations.  However,  it  is  certain  that  the  outgassing 
flux  does  decrease  over  time,  and  this  is  not  observed  in 
the  data.  Additionally,  it  is  hard  to  conceive  that  the 
satellite  potenials  are  such  that  the  randomly  ionized  flux 
is  focussed  into  a  beam  one  energy  channel  wide  at  the 
detector.  Moreover,  it  is  unreasonable  that  the  focussing 
should  be  such  that  a  small  change  in  the  detector  look 
angle  yields  such  a  definite  change  in  the  observed  energy, 
as  in  the  phenomena  of  the  triangle  peaks.  In  sum  then,  it 
does  not  appear  likely  that  ionization  of  outgas  products  is 
the  source  of  the  observed  flux. 


When  this  work  began,  it  was  initially  assumed  that 
the  ion  source  was  outgassing.  As  the  work  progressed,  it 
became  apparent  that  there  were  serious  drawbacks  in  this 


assumption,  as  noted  above 


While  searching  for 


alternatives,  the  advantages  of  sputtering  as  a  source 
became  clear.  First,  besides  a  small  dose  dependence,  there 
is  no  reason  for  the  sputtering  yield  to  decay  over  time. 
This  attribute  was  very  important  in  light  of  the 
significant  shadow  peaks  observed  in  1981.  Second,  the 
emission  of  sputtered  particles  from  a  differentially 
charged  surface  would  lead  to  observed  beams  at  the 


detectors . 


As  the  particles  left  the  surface  with  their 


various  initial  angles  and  energies,  they  would  be 
accelerated  by  the  existing  electric  field  into  an  almost 
parabolic  orbit.  Thus,  for  each  energy,  only  particles 
following  specific  trajectories  would  reach  the  detector, 
limiting  observation  angles.  Since  the  detector  energy 
channels  have  a  resolution  of  20  %,  the  bandwidth  of  the 
channels  increases  with  energy.  Therefore,  a  peak  which 
would  appear  broad  when  distributed  over  the  lower  energy 
channels,  can  be  compressed  into  a  narrow  peak  in  the  high 
channels.  The  low  and  high  energy  channels  of  the  SCATHA 
detectors  are  listed  in  Table  8.  The  similarity  between  the 
observed  data  and  the  data  from  the  Minnesota  spots  also 
became  evident.  It  seemed  reasonable  that  if  such  a  process 
could  occur  for  secondary  electrons,  it  could  occur  for  what 
were  essentially  "secondary"  ions.  Additionally,  the 
sputtering  theory  could  explain  the  diffuse  spectrum  as  the 
receipt  of  particles  sputtered  from  a  non  differentially 
charged  surface.  Without  the  charging,  the  ions  are  not 
focussed  into  the  higher  channels,  and  can  appear  as  a  broad 
low  energy  distribution.  The  third  phenomena,  the  shadow 
peak,  can  be  described  by  two  mechanisms.  The  first  is  a 
motionless  detector,  not  scanning  in  angle,  and  observing 
one  trajectory,  i.e.,  one  energy.  Alternatively,  the 
detector  could  be  scanning  laterally  across  an  emitted  beam. 
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as  opposed  to  vertically,  and  thus  be  picking  up  the  same 


energy  at  each  angle.  Figure  30  illustrates  the  geometry. 

The  problem  with  sputtering  as  a  source,  was  that  it 
had  never  seemed  to  be  able  to  provide  an  adequate  flux  of 
ions.  However,  the  only  previous  mention  of  sputtering  had 
been  prior  to  the  advent  of  the  Sigmund  theory,  and  before 
it  was  acknowledged  that  the  magnetospher ic  regions 
contained  large  numbers  of  oxygen  ions  {Ref.l}.  Preliminary 
calculations  at  this  point  indicated  that  the  flux  due  to 
oxygen  sputtering  was  sufficient  to  account  for  the  observed 
phenomena.  It  then  remained  to  model  the  situation  and 
determine  more  accurately  the  flux  at  the  detectors  due  to 
sputter ing . 

It  was  now  naressary  to  ev’aluate  the  physical 
characteristics  of  the  environmental  interaction  with  the 
spacecraft.  As  the  bulk  of  the  vehicle  surface  was  covered 
with  solar  panels,  the  most  likely  target  of  the  incident 
ions  would  be  amorphous  silicon  dioxide.  For  a  magnetically 
disturbed  time  period,  the  ion  flux  would  be  composed 
primarily  of  oxygen  and  hydrogen.  The  Sigmund  theory  has  a 
strong  dependence  on  the  mass  of  the  incident  ions,  and  it 
is  reasonable  to  expect  a  much  larger  sputtering  yield  for 
the  heavier  ion.  See,  for  example.  Figure  14  for  and  0"^ 
on  amorphous  silicon. 

The  next  step  was  to  use  the  equations  (1-7)  and  (1- 
8)  to  calculate  the  yield.  One  roadblock  to  this  was  the 


Detector 

Scamma  vertically  TTirough  Energy 


Detector  Scanrttig  Across  Angle  Receiving 
Constant  Energy 


inability  of  the  Sigmund  theory  to  adequately  treat 
sputtering  from  multi-atomic  targets.  There  are  some  data 
points  available  for  impacting  molecular  oxygen  on  silicon 
and  silicon  dioxide,  and  they  are  listed  in  Table  9 
{Ref. 36}.  If  the  data  from  Table  7  are  applied  to  Table  9, 
we  see  that  the  sputtering  yield  for  0*  on  silicon  will 
range  from  .1  to  .18  at  10  keV.  Similarly,  the  yield  on 
silicon  dioxide  will  range  from  .23  to  .44  atoms  per  ion  at 
10  keV.  The  sputtering  yield  calculated  from  the  Sigmund 
theory  for  atomic  oxygen  on  silicon  is  too  high,  and  more 
accurately  approximates  the  yield  from  SiOz.  We  will, 
however,  take  the  worst  case  approach,  and  calculate  the 
yield  for  silicon.  To  do  this,  we  will  increase  the  binding 
energy,  to  give  us  values  more  in  line  with  those  from  the 
experiment.  The  two  yield  curves  are  displayed  in  Figure 
31.  The  new  binding  energy  is  26.0  eV. 

Given  that  this  new  curve  now  reasonably  describes 
the  sputtering  yield,  the  problem  is  to  calculate  how  many 
particles  emitted  from  a  surface  will  reach  a  differentially 
charged  detector  for  a  realistic  incident  flux.  The  total 
yield  of  sputtered  particles  will  be  described  by 

(3-1)  ^  ^  I"" 

In  this  expression,  r  is  the  maximum  energy  a  sputtered 
particle  can  have,  that  is,  the  maximum  amount  of  energy 


transferred. 


This  maximum 


13 


calculated  by 


Sputtering  Yield  vs.  Energy,  0+  on  Si 


and  is 


r  =  4  M:l  M2  Em, 

(Mi+Ma ) = 


equivalent  to  the  maximum  energy  transferred  in  a 
binary  collision.  The  function  f(v)  is  the  distribution 
function  of  the  incident  ions,  and  S(v,9)  is  the  sputtering 
yield . 

With  the  total  yield  in  hand,  we  may  now  attempt  to 
calculate  the  countrate  at  the  detector.  The  total  number 
of  counts  the  detector  receives  may  be  described  by  the 
following  equation. 


(3-2) 


Counts  =  fdt  Pd  A  Pf(v)  (^‘'n)d^v 


This  assumes  that  the  detector  is  differential  in  time, 
area,  solid  angle,  and  energy.  If  we  integrate  over  the 
area  increment  and  time,  we  get. 


■3) 


Counts=T  3A 


r”r 


'f  (v)  ( V'-ft)dv 


Thus,  the  countrate  is 


(3-4) 


Countrate=Counts/Time=(5A 


JJ 


v*f  ( v)  (^»  n)dvdil 


To  continue,  we  must  assume  that  limits  of  integration  on  dft 
are  not  a  function  of  energy,  and  that  the  value  of  the 
integral  may  be  given  by  a  simple  number,  characteristic  of 
the  detector.  This  is  not  actually  true,  as  the  limits  on 
integration  depend  strongly  on  energy,  since  f(v)  may  not 
fill  the  entire  5  by  7  degree  aperture  of  the  detector.  The 
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inaccuracy  of  this  assumption  will  surface  later, 
assumption,  we  derive. 


With  the 


(3-5)  Countrate  =  CR  =  (3  A  <3ft  Jv=  f  ( v )  ( v- n  )  dv 

The  two  constants,  (3A  and  6Si,  are  combined  into  the 
geometric  factor  G,  a  constant  for  the  detector.  For  our 
detectors,  G=1.6  *  m^  ster.  Due  to  the  differential 
aperture  of  the  detector,  the  value  of  the  dot  product  n 
is  approximately  v  (i.e.  cos  Dal). 

We  must  now  determine  the  distribution  function  of 
the  emitted  particles.  From  the  earlier  discussion  of  the 
Thompson  Theory,  it  is  evident  that  the  energy  distribution 
of  the  yield  flux  is  proportional  to  E/CE+Uo)^.  That  is,  we 
have  the  function 


( 3-6 )  dY  =  A  E 

dE  (E+U0)=> 

The  flux  Y,  is  also  given  by 
(3-7)  Y=BJ“Ef(E)dE 

where  f(E)  is  the  distribution  function  as  a  function  of 
energy  of  the  emitted  particles,  and  B  is  a  constant  that 
involves  the  integral  over  the  solid  angle.  If  we  combine 
these  two  relationships,  it  is  evident  (disregarding 
constants)  that 

(3-8)  f(E)  =  I  dY  =  1 

E  dE  (E  +  Us, ) 
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We  will  now  determine  the  appropriate  normalization 


constants.  From  the  definition  of  the  distribution 
function. 


(3-9  ) 


f  (  V )  d  ^  V 


Continuing  the  calculation,  we  assume  that  the  yield  flux 
and  the  distribution  function  have  similar  angular 


dependencies.  Thus,  f ( v, 0 ) =f ( v ) cos9 .  If  we  break  down  d^v, 
and  substitute  for  energy,  we  get 


( 3-10) 


n=AJd0js i 


sine  cose  d9Pf  (E)  f  2E  ) '^ 


In  the  case  of  emission  from  a  surface,  #  will  range  from  0 
to  2k,  and  0  will  vary  from  0  to  k/2.  Integrating  over 
these  limits,  the  equation  reduces  to 


(3-11)  n=Ait(2/m=»)'* 


r _ e:i_ 

J  (E  +  Us 


This  integral  will  extend  the  full  energy  range  of  the 
emitted  particles,  that  is,  from  0  energy  to  r,  the  maximum 


energy  transferred  in  the  collisions. 


If  one  then 


substitutes  a  dummy  variable,  x=,  for  E  the  integral  has  a 
standard  form  {Ref.  37}.  Solving,  we  then  get 


n=Am(2/m'’)'*  P*" 


r _ s: 

J0  (E  + 


n=ATC(8/m=’)’*  P- 


P*  x^ _  c 

Je  (X^  +  Ub)^ 


s  = 


n  m  i  « 


n=Ait(8/m=)''  {(-x)/{4(x*  +  U® )  =  )  +  x/(8Ub(x=  +  Ub)^) 
+  (1/8(Ub'>)'*)  tan-^  (x/CU^)'*)  }-8 

n=ATt(8/m=’)'*  {(r'*/(8(U0  +  r)))(l/U0  -  2/{U<a  +  D) 

+  {1/(8U»=’^*)  )  tan-=>-  (r/UB)"} 

Defining  8  as  the  large  constant  on  the  right  side,  ve  get 

n=ATt  (  8/in^  ) '*  8 
or 

(3-12)  A=n  m3^  =  /Tr88'* 

Using  this  normalization  constant,  the  distribution  function 
becomes 

(3-13)  f(E,e)=  (n  m^-^  =  /Tt08'‘)  cos0  1 

(E  +  U»)= 

If  we  then  calculate  the  yield  flux,  Y,  we  obtain 

Y=(l/2m)’*  (n/0)  E  dE 

U  (E  +  Ub)* 

(3-14)  Y=(l/2m)'‘  (n/0){-(r  +  U,/2)/(r  +  Ub  )  =  +  (1/2Ub)} 

Y=(l/2m)'‘  (n/0)  <3 
or 

(3-15)  n=(  2m)’*0Y/(3 

Therefore,  substituting  in  the  distribution  function, 

f(E,0)  =111  1  COS0 

2n  <5  (E  +  Ub)3 


(  3-16  ) 


We  now  have  an  expression  for  the  distribution 
function  of  the  emitted  particles.  This  expression  can  then 
be  inserted  with  one  modification  into  eqn .  (3-5)  to 

determine  the  countrate  at  the  detectors.  The  necessary 
modification  accounts  for  the  differential  charge  between 
the  emitting  surface  and  the  detectors.  For  the  particles 
to  be  accelerated  to  the  detector,  the  surface  must  be 

charged  positively  relative  to  the  detectors,  therefore  « 
will  be  negative.  Since  E  will  map  to  E  +  (Liouville’s 

Theorem),  we  will  obtain 

(3-17)  f(E,0,«)=  itil  Y  _ 1 _  COS0 

2Tt  (3  (  (E-«)  +  Ub)  ’ 

Here,  we  have  assumed  that  the  ions  are  singly  charged,  and 
I#  will  be  converted  to  eV  from  V  on  this  basis.  The 
countrate  expression  is  then, 

(3-18)  CR=G(2/m*)  J  f(E,0,«)  E  dE 

(3-19)  CR  =  2GY  COS0  p  _ E _  dE 

Tt<5  J  (  (E-«)  +  Ub  )  = 

To  find  the  countrate,  we  must  now  integrate  for  each  of 

the  detector  energy  channels  over  its  energy  range.  As  this 

resolution  is  20%,  the  upper  and  lower  limits  are  1 . lEc  and 
0.9Ec,  where  Ec  is  the  central  energy  of  the  energy  channel. 
The  integral  in  eqn.  (3-19)  is  a  standard  form, 
countrate  may  finally  be  written  as. 


and  the 


then  the 
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If  the  energy  channel  limits  are  both  less  than  0, 
countrate  over  that  channel  is  zero.  If  the  lower  limit  is 
below  0,  then  the  integral  is  taken  between  the  upper  limit 
and  0 . 


B.  CALCULATION  OF  THE  SPUTTERING  FLUX  AT  THE  DETECTOR 

Using  the  Sigmund  Theory  developed  above,  and  the 
derived  expression  for  the  countrate  at  the  detector,  we 
developed  a  program  to  calculate  the  sputtering  yield, 
incident  and  emitted  flux,  and  the  reponse  of  the  detector 
to  changes  in  both  the  energy  channel  and  differential 
charge,  0.  Data  for  the  ambient  plasma  compositions  was 
available  from  the  Lockheed  Ion  Mass  Spectrometer  for  days 
83,  86,  and  200  of  1979.  The  values  utilized  are  listed  in 
Table  10. 

The  first  portion  of  the  program  was  the  calculation  of 
the  total  yield  or  emitted  flux.  This  value  is  calculated 
from  the  following  expression; 


3-21)  Y  =  P-  V  f(v)  S(v,0)  d’v 


In  this  relation,  v  f(v)  is  the  flux  of  the  ambient  ions  as 
a  function  of  velocity,  and  S(v,0)  is  the  sputtering  yield 
as  a  function  of  velocity  and  angle.  The  factor  S„(£)  in 
S(v,0)  exists  in  tabulated  form,  and  is  not  described  by  a 
simple  continuous  function.  If  the  points  are  fitted  to  a 


curve,  we  arrive  at  the  equation; 

(3-22)  S„(e)=  .35874  -  (. 17885 ) loge  - (  . 1 4 3 59  )  ( loge  )  = 

+( .06193)(log£)^  f  (. 04228 )( logs ) * 

-( .00541) {log£)=  -( .00630) (log£)“ 

-  I  .  00092  )  ( log£  ) 

This  equation  and  the  available  tabulated  points  are 

displayed  in  Figure  32.  Since  Sn(£)  is  a  complicated 
function,  the  integral  was  evaluated  numerically.  The 

function  was  calculated  at  a  thousand  points  over  an  energy 
range  that  depended  on  the  temperature  of  the  ambient 
Maxwellian  distribution.  For  a  4000  eV  distribution,  the 
high  energy  was  50000  eV  and  the  low  energy  was  100  eV.  100 
eV  was  chosen  as  the  lower  limit,  as  most  sputtering  yields 
are  approximately  zero  by  that  point.  At  the  50000  eV 
point,  the  value  of  the  distribution  function  has  reduced  by 
a  factor  of  10“^"  from  the  peak  value. 

The  yield  obtained  in  this  method  was  then  multiplied  by 
.02,  as  approximately  2%-4%  of  the  emitted  flux  is  ionized 
{Ref.  29}.  This  ion  yield  is  substituted  in  the  countrate 
integral  to  obtain  the  detector  responses.  The  values  of 
the  integrals  were  varied  over  energy  channel  and 

differential  charge  to  determine  if  the  responses  were  very 

narrow  as  observed.  Figures  33  shows  the  values  of  the 
incident,  emitted,  and  ion  fluxes  for  oxygen  as  a  function 
of  energy  for  day  86  of  1979.  The  figures  for  days  83  and 
200  are  similar.  From  this  it  is  very  clear  that  the 
incident  flux  is  much  greater  than  the  emitted  ion  flux. 
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Thus  there  is  a  net  accumulation  of  positive  charge  on  the 
spacecraft  due  to  the  ambient  ion  current. 

Figures  34  and  35  are  displays  of  the  countrate  as  a 

function  of  energy  channel  for  a  constant  (35  eV) 

differential  charge,  for  days  86  and  200  for  oxygen  induced 
emission.  We  see  that  the  countrate  is  high  in  3-4  of  the 
energy  channels.  The  data  indicates  that  it  in  fact  should 
be  concentrated  in  only  one  channel.  However,  the  above 
equations  assume  that  the  detector  surface  is  a  flat  plate, 
and  that  all  particles  reaching  that  plate  are  counted. 
This  is  in  fact  not  true.  The  detector  viewing  cone  is 
probably  not  completely  filled,  as  noted  previously,  and  we 
are  seeing  the  error  due  to  treating  the  factor  6Si  as  a 

constant  in  eqn.  (3-3).  This  factor  will  serve  to  reduce 

the  countrate  overall,  and  since  particles  at  the  edges  of 
the  viewing  cone  (with  different  energies  from  the  cone 
center)  will  be  reduced,  the  peaJ<  will  most  likely  become 
narrower.  The  response  for  a  specific  channel  and  varying 
differential  charges  is  quite  similar.  That  is,  a  given 
channel  will  only  respond  for  a  certain  narrow  range  of 
differential  charge.  Similar  calculations  for  the  ambient 
hydrogen  population  on  the  same  days  indicate  that  it  can 
also  generate  a  significant  flux  at  the  detectors,  simply 
due  to  the  large  amount  present.  The  relative 
contributions  of  these  two  ambient  plasma  components  very 
widely,  depending  on  the  temperature  and  density  of  their 


respect i ve 


d istr ibut ions . 


general , 


the 


oxygen 


contribution  is  dominant. 

The  countrates  from  the  detectors  are  available  for  the 
above  three  days.  On  day  83  of  1981,  the  measured 
countrates  varied  from  29  -  800  counts  per  second,  with  the 
average  at  360.  For  day  86,  the  range  was  24  -  1080,  with 
the  average  at  300.  Day  200  was  much  lower,  between  25  and 
213,  with  an  average  value  of  120.  The  available  mass 
spectrometer  data  was  averaged  over  one  hour  periods,  so  we 
will  generally  calculate  only  one  representative  number  for 
the  charging  time  frame.  The  calculated  countrates  were 
380,  325,  and  400  for  days  83,  86,  and  200,  respectively. 
Considering  the  approximate  nature  of  the  calculations,  the 
theoretical  responses  compare  favorably  with  the  measured 
data,  with  the  exception  of  day  200.  However,  it  has 


already  been  noted  that  this  data  from  the  Hi  detector 
indicated  a  sputtering  source  different  from  those  seen  on 


the  Lo  detector.  This  further  supports  that  observation,  as 
the  sputtering  yield  for  this  event  is  clearly  much  lower. 


As  the  assumptions  in  this  calculation  tend  to  minimize 
the  number  of  detected  particles,  it  seems  reasonable  to 
conclude  that  sputtering  from  a  differentially  charged 
surface  is  sufficient  to  account  for  the  narrow  triangle 


peaks.  These  discrete  responses  will  serve  to  explain  the 
triangle  and  shadow  peaks.  This  leaves  the  diffuse 
background.  We  require  a  broad  response,  with  the  particles 


scattered  equally  across  the  energy  channels.  Figure  36 
gives  the  log  distribution  function  as  a  function  of  energy 
for  an  emission  angle  of  0  degrees.  The  distribution 
function  has  its  maximum  in  the  low  energy  channels,  and 
numerically  compares  favorably  with  Figures  21  and  24.  The 
broad  spectrum  may  then  be  described  as  the  returning 
emissions  from  non-d i f f erent ial ly  charged  surfaces. 

We  must  now  explain  the  shadow  peaks.  These  appear  to  be 
a  combination  of  two  occurrences,  acceleration  from  a 
differentially  charged  surface  to  a  motionless  detector,  or 
acceleration  from  a  differentially  charged  surface  to  a 
detector  scanning  laterally  across  the  emitted  beam,  as 
opposed  to  vertically.  See  Figure  31  for  a  view  of  the 
geometry.  Since  this  is  then  a  variation  of  the  triangle 
peaks,  the  source  and  results  would  be  the  same. 

The  trajectories  to  the  detectors  can  be  recreated  by 
the  NASCAP  program.  Figure  37  displays  the  particle 
trajectories  for  the  observed  angles  and  energies  at  the  Lo 
detector  on  day  86.  From  top  to  bottom,  the  lines  are  as 
follows;  a)84  eV,  145  degrees,  or  12.5  degrees  above  the 
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Figure  36.  Channel  Response  for  a  Zero  V  Differential 
Charge 
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the  sputtered  particles  follow  specific  trajectories  in  the 
electric  fields.  The  limited  view  cone  of  the  detector  only 
picks  up  specific  energies. 

NASCAP  was  also  run  for  the  potentials  observed  on  day 
94.  The  code  was  run  to  simulate  the  responses  of  the  Lo, 
Hi,  and  Fix  detectors  to  a  range  of  incoming  particle 
energies.  The  energies  ranged  from  2.0  to  90.1  eV.  As 
noted  above  in  the  observations,  the  spacecraft  was  charged 
to  -60  V,  and  peaks  were  noted  at  various  lower  energies  on 
the  detectors.  The  Hi  detector  was  parked  at  67  degrees, 
and  the  Lo  at  91  degrees.  Thus,  Hi  and  Lo  detectors  were 
sampling  approximately  the  same  region  of  space. 

Figures  38  and  39  depict  the  Lo  detector  trajectories  in 
two  different  planes.  The  first  shows  the  top  of  the 
spacecraft  surface.  Particles  at  energies  less  than  43  eV 
are  returned  to  the  spacecraft.  The  trajectories 
corresponding  to  observed  particles  are  marked  with  an 
asterisk.  The  second  figure  is  a  complementing  view  from 
the  left  side  of  the  first  figure.  These  two  figures 
indicate  that  the  actually  observed  particles  are  emitted 
from  the  satellite's  forward  end.  Additionally,  there  are 
gaps  from  which  no  particles  are  seen  at  the  detector.  An 
explanation  for  this  is  the  existence  of  patches  of 
different  nonconducting  material  on  the  forward  surface  for 
use  in  the  Satellite  Surface  Potential  Monitor  (SSPM).  When 
these  patches  (kapton,  silvered  teflon,  and  quartz)  charge 


differentially,  the  particles  reaching  the  detector  from 
each  will  follow  different  trajectories,  and  have  different 
energies.  Additionally,  the  S-band  antenna  mast  on  the 
forward  end  is  nonconducting,  and  may  be  a  particle  source. 

Figures  40  and  41  are  similar  views  for  a  model  of  the 
Ki  detector  trajectories.  We  see  that  the  origin  of  the 
particles  is  also  on  the  upper  surface,  and  the  data  shows 
the  same  striations.  Again,  the  only  nonconducting  surfaces 
on  the  forward  end  are  the  SSPM  and  the  antenna.  The  detail 
on  NASCAP  is  not  sufficient  to  study  the  extremely  strong 
local  electric  fields  expected  in  this  region.  If  the 
satellite  potential  is  100  V,  then  the  field  due  to  the 
sheath  is  approximately  1  V/m.  In  the  forward  surface 
vicinity,  with  a  differential  potential  of  30  V,  the  fields 
will  be  on  the  order  of  30-60  V/m,  and  will  significantly 
distort  the  local  trajectories  of  particles.  Thus,  the 
NASCAP  trajectory  plotting  routine  is  somewhat  unreliable  in 
this  region.  Looking  at  both  figures,  however,  the  origins 
for  the  two  detectors  are  reasonably  close,  and  allowing  for 
the  local  fields,  may  possibly  come  from  the  same  surface. 
As  a  working  hypothesis,  we  will  assume  that  the  SSPM  is  the 
source  of  the  emitted  particles. 

Figures  42  and  43  model  the  fixed  detector.  The  views 
are  similar  to  those  above.  We  see  that  the  particles 
viewed  by  the  fix  detector  come  from  completely  different 
regions  of  the  spacecraft.  All  the  obviously  reattracted 


97 


12345578  9  10  11  1  2  1  3  1  4  15  1  6  17 


T-AXIS 


Figure  43.  N\SCAP  Trajectory,  Day  94,  Fix  Detector 
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particles  have  energies  of  less  than  43  eV,  as  did  the  other 
two  detectors.  We  can  see  that  the  trajectories  of  the 
particles  skip  over  the  bellyband  regiou.  This  is 
reasonable,  as  the  bellyband  would  be  expected  to  be  at  the 
satellite  mainframe  potential,  and  thus  particles  emitted 
would  not  be  accelerated  by  a  differential  charge.  The 
trajectory  calculations  for  the  Fix  detector  also  support 
the  theory  of  sputtering  from  differentially  charged 


It  is  reasonable  to  conclude  that  sputtering  is  a  likely 


source  of  the  observed  phenomena.  It  can  satisfactorily 
explain  the  flux  and  structure  of  the  triangle  and  shadow 
peaks.  It  also  can  account  for  the  observed  broad  spectrum 
of  ions.  These  fluxes  are  large  enough  to  be  measured  by 
the  detectors,  but  are  they  large  enough  to  cause 
substantial  damage  to  the  spacecraft  surface? 

For  this  data  set,  the  sputtering  yield  was  less  than 
the  incident  flux,  so  that  the  net  result  was  particle 
accretion  on  the  surface.  Coincidentally,  the  data  in  this 
thesis  were  taken  from  relatively  low  energy  days.  On 
higher  energy  days,  the  sputtering  yield  increases,  and  the 
net  result  is  surface  depletion.  However,  the  surface 
damage  caused  by  this  process  appears  to  be  minor.  Assuming 
a  higher  energy  day,  with  emission  on  the  order  of  twice  the 
incident  flux,  one  atomic  layer  will  be  removed  from  a 
silicon  dioxide  surface  every  8  years.  Since  the  solar 
cells  are  covered  by  the  glass  shields,  they  will  be 
impervious  to  the  surface  sputtering. 

Another  facet  of  this  is  the  damage  to  surface 
conducting  coatings  on  spacecraft.  These  coatings  are  only 
several  atoms  thick  and  may  be  expected  to  sputter  from  the 
surface  within  several  months  of  the  satellite  launch.  This 
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may  account  for  the  observed  charging  on  the 
satellite. 
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Since  the  energy  spectrum  of  the  sputtered  particles 
dies  off  quite  rapidly,  and  the  observed  differential 
charges  are  usually  not  large,  it  is  not  likely  that  any 
great  contami nat i on  of  the  charging  peak  results.  The 
sputtering  in  general  will  probably  be  uniform  over  a  given 
material.  Surfaces  will  be  worn  down  over  large  areas  as 
opposed  to  experiencing  significant  localized  damage. 
Therefore,  it  appears  that  no  short  term  damage  will  result 
from  sputtering  on  spacecraft. 

Further  work  is  required  to  accurately  calculate  the 
sputtering  yields  from  various  surfaces  to  prove  or  disprove 
these  conclusions.  At  this  point,  given  the  current  state 
of  the  sputtering  field,  it  is  not  possible  to  state  with 
any  certainty  whether  extensive  damage  will  result  from 
ambient  ion  bombardment.  There  appears  to  be  no  work  in 
progress  anywhere  to  determine  experimentally  the  sputtering 
yields  for  most  ion-target  combinations.  Since  the  Sigmund 
Theory  is  so  limited,  these  results  are  necessary  for  any 
actual  calculations  of  the  effect  of  the  ambient  heavy  ion 
bombardment  on  space  vehicles. 

This  thesis  demonstrates  that  sputtering  is  the  probable 
cause  of  the  observed  low  energy  ion  fluxes.  Additional 
work  is  required  to  refine  the  NASCAP  trajectory  information 
and  couple  it  with  the  theory  developed  here  to  model  the 
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STORM  TIME  RING  CURRENT  COMPOSITION  {Ref.  8} 


Date 

07  July 

26  Nov 

09  Mar 

3  0  Me 

1977 

1977 

1978 

1978 

H-^  Flux 

10'’cm~  =  s~^ 

1.26 

4.20 

13 . 0 

6 . 0 

Dens  i  ty 

10-=cm-3 

H-' 

11 

44 

143 

58 

He* 

.72 

4  .  3 

10.9 

6.6 

0* 

8.3 

91 

228 

107 

Mean  Energy 
keV 

H* 

8.3 

5.7 

5.1 

6.4 

He* 

4.5 

4 . 3 

5.4 

2.7 

0* 

4.9 

3.8 

5.3 

4  .  4 

% 


44^ 
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STORM  TIME  MAGNETOSPHERE  COMPOSITION  {Ref.  8} 


Date 

04  May 
1977 

29  July 
1977 

29  Jan 
1978 

02 

197 

Flux 

10‘'cm“*s 

3 . 3 

-  a. 

2.4 

10 

2 . 8 

Dens i ty 

26 

20 

10 

27 

He-^-^ 

1.5 

.  16 

9 . 3 

.  37 

He’" 

1.8 

.9 

.  3 

.79 

O’" 

3.0 

13 

4.9 

60 

Mean  Energy 
keV 

H-" 

9.6 

8.5 

4.2 

6.8 

He'"’" 

20.3 

15.9 

8.6 

14 

He-" 

4.8 

3.9 

6.4 

4 . 8 

0" 

4.5 

4.7 

5.4 

4.1 
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VALUES  OF  THE  REDUCED  STOPPING  POWER  C Ref. 27} 


CALCULATED  AND  MEASURED  SPUTTERING  YIELDS  FOR  VARIOUS 
MATERIALS  {Ref.  23} 

THE  NUMBERS  IN  PARENTHESES  ARE  THE  YIELDS  CALCULATED 
FROM  THE  SIGMUND  THEORY 


TARGET 

Soutterina  Ratio  (atoms/ion) 

MATERIAL 

UaCeV) 

Zt 

Ne- 

Ar-^ 

Kr:-^ 

Xe* 

Pb 

2.01 

82 

3.6 

10 . 5 
(29.7) 

24.0 

(44.0) 

44 . 5 
(74.6) 

Ag 

2.94 

47 

4 . 5 

10.8 

(12.7) 

23 . 5 
(20.1) 

36 . 2 
(24.0) 

Sn 

3.11 

50 

1.8 

4  .  3 

(12.5) 

8 . 5 
(19.5) 

11.8 
(  4.9) 

Cu 

3.46 

29 

3.2 

(3.7) 

6.8 

(6.7) 

11.8 

(11.9) 

19.0 

(15.5) 

Au 

3.79 

79 

3.6 

10,2 

(16.1) 

24 , 5 
(23.9  ) 

39 . 0 
(27.8) 

Pd 

3 . 87 

46 

2.5 

5.3 

(9.24) 

10.5 

(14.1) 

14 . 4 

(18.1) 

Fe 

4 . 29 

26 

1.3 

2.3 

(5.15) 

4.0 

(8.92) 

4.9 

(11.7) 

Ni 

4.43 

28 

1.4 

3.5 

(5.35) 

5.6 

(9.01) 

7.6 

(11.8) 

V 

5.3 

23 

0.3 

1.0 

(3.85) 

1.7 

(6.21) 

1.9 

(8.9) 

Pt 

5.82 

78 

1.9 

5.3 

(9.35) 

11.3 

(15.3) 

16 . 0 
(19.2) 

Mo 

6 . 82 

42 

0.6 

1.5 

(5.14) 

2 . 7 
(7.70) 

3.8 

(10.1) 

Ta 

8.06 

73 

0.7 

1.6 

(6.62) 

3.1 

(10.0) 

4.0 

(12.3) 

W 

8.70 

74 

1.0 

2 . 3 

(5.07) 

4.7 

(9.31) 

6 . 4 

(11.9) 

TABLE  $ 
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RATIO  OF 
AND 

SPUTTERING 
ATOMIC  ION 

YIELD  PER  ATOM  FOR 
BOMBARDMENT  {REF. 

MOLECULAR 

31} 

Project!  les 

Si 

Ag 

Au 

C1-C1  = 

•  •  • 

1 . 09 

.  .  . 

Se-Sea 

1.15 

1.44 

1.44 

Te-Tez 

1 . 30 

1 .67 

2 . 15 

table  8 


SCATHA  PARTICLE  DETECTOR  ENERGY  CHANNELS  {Ref.  4} 


Lo  Detector 


Step 

Energy ( eV ) 

Step 

Energy ( eV ) 

0 

-0.334 

32 

26 . 26 

1 

-0.29 

33 

30.10 

2 

-0.23 

34 

34 . 57 

3 

-0.16 

35 

39 . 70 

4 

-0 .08 

36 

45.55 

5 

0.00 

37 

52 . 22 

6 

0.09 

38 

60.00 

7 

0 . 22 

39 

68 . 77 

8 

0.34 

40 

78.55 

9 

0.50 

41 

90.10 

10 

0 .67 

42 

103.32 

11 

0.87 

43 

118 .32 

12 

1.10 

44 

135.65 

13 

1.36 

45 

155.43 

14 

1.66 

46 

178.09 

15 

2.00 

47 

203.98 

16 

2.39 

48 

233.42 

17 

2.83 

49 

267 . 20 

18 

3.34 

50 

306 . 30 

19 

3.94 

51 

350.74 

20 

4.60 

52 

401.96 

21 

5.37 

53 

457.73 

22 

6.25 

54 

525.50 

23 

7.25 

55 

601.05 

24 

3.40 

56 

687.71 

25 

9 . 72 

57 

787.70 

26 

11.23 

58 

902 . 13 

27 

12.96 

59 

1033.23 

28 

14.96 

60 

1184.33 

29 

17.23 

61 

1355.42 

30 

19.83 

62 

1529 .85 

31 

22.80 

63 

1777.60 

TABLE  8 


Wil' 


n.in 


CONTINUED 


Hi  Detector 


Step 

Energy  (eV) 

Step 

Energy 

0 

-4.68 

32 

1205.49 

1 

-2 .35 

33 

1383  .  30 

2 

0.32 

34 

1586.90 

3 

3  .  38 

35 

1820 . 02 

4 

6 .88 

36 

2086.94 

5 

10 . 89 

37 

2392 . 57 

6 

15.49 

38 

2742 . 51 

7 

20.75 

39 

3143 . 19 

8 

26 .77 

40 

3601.97 

9 

33.66 

41 

4127 . 28 

10 

41.56 

42 

4728.76 

11 

50 .60 

43 

5417 . 45 

12 

60.95 

44 

6206.00 

13 

72.80 

45 

7108 . 89 

14 

86 .37 

46 

8142.70 

15 

101.91 

47 

9326 .41 

16 

119.71 

48 

10681.76 

17 

140.08 

49 

12233.63 

13 

163.41 

50 

14010.53 

19 

190.12 

51 

16045 . 08 

20 

220.70 

52 

18374.64 

21 

255.72 

53 

21041.98 

22 

295.81 

54 

24096 . 09 

23 

341.72 

55 

27593.05 

24 

394.29 

56 

31597.06 

25 

454.45 

57 

36181.66 

26 

523.39 

58 

41431.02 

27 

602.30 

59 

47441.54 

28 

692.65 

60 

54323 . 58 

29 

796.11 

61 

62203 . 53 

30 

914.56 

62 

71226.06 

31 

1050.19 

63 

81556 . 86 

TABLE  Ig 


V' 

JUS 

'•ft 

rii 

•'ti'j 


AMBIENT  PLASMA  VALUES  FOR  DAYS  83,  86,  AND  200  OF  1979. 
DATA  FROM  THE  LOCKHEED  ION  MASS  SPECTROMETER. 


Particle 


Number  Density  (m' 


Energy ( keV ) 


83/ 

H-t- 

3 .845 

* 

10  = 

8 .852 

1979 

H=* 

2.487 

* 

10  = 

120.778 

He 

2.865 

* 

10  = 

33.992 

0-" 

3.188 

* 

10  = 

4  .  217 

86/ 

H-^ 

2.606 

* 

10  = 

5 .889 

1979 

Ha-" 

2.850 

* 

10  = 

43.728 

He 

2 . 410 

* 

10  = 

3 .662 

0- 

1.635 

* 

10  = 

4 .035 

200/ 

H- 

9.064 

* 

10  = 

5.271 

1979 

4.209 

* 

10'‘ 

7 . 279 

He 

2.169 

* 

10  = 

0.420 

0-^ 

4.093 

10  = 

1.957 
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